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SECTION 1 . 


INTRODUCTION AND SOMMAR!^ 


1.1 INTRODUCTION 

This report has been prepared as the Final Report in a series 
of studies performed by the Charles Stark Draper Laboratory, Inc. , 

(CSDL) , for the NASA Goddard Space Flight Center (GSFC) , during the 
period from 13 April 1973 to 31 December 1975. The present report 
gives the results of a recent study conducted during the period from 
1 February 1974 to 31 December 1975 which is referred to as the Arti- 
ficial Landmark Implementation Study (ALIS) . This study was an inves- 
tigation into the feasibility, practicality, and cost of establishing 
a national system of artificial landmarks suitable for automated recog- 
nition in the electronic processing of imagery data from an Earth 
Observation Satellite (EOS) . The results of the other studies performed 
rnider this contract are presented in References 1, 14, 15 and 16. 

1.1.1 Background 

The ALIS study was motivated by the preliminary results of 
another study performed by CSDL for NASA (GSFC) , during the period from' 

1 February 1974 to 31 July 1974. That study, herein referred to as the 
Landmark Utilization Study (LUS) , was an investigation into the feasibil- 
ity of using ground truth data, associated with selected landmarks in 
the EOS imagery, and also star measurement data in combination with 
landmark data, to estimate the orbital ephemeris, attitude and gyro bias 
drift of EOS satellites. 

The preliminary findings of the LUS study indicated that the 
automated recognition of a small number of ground targets (in the EOS 
imagery data) in combination with star measurement data (from an on- 
board star tracker) could be used to determine the orbital ephemeris of 
the EOS spacecraft and to create a significant improvement in the geo- 
metric mapping accuracy of the telemetered imagery data. These prelimi- 
nary findings of the LUS study were subsequently substantiated by a more 
detailed analysis which indicated that the above benefits could be ob- 
tained with as little as ten star sightings per orbit and just one or 
two landmark sightings per orbit. These findings were derived by 
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means of a statistical covariance analysis, employing optimal filtering 
and smoothing techniques to process the star and landmark measurement 
data along with the ordinary attitude and gyro bias drift data. In this 
way, star and landmark data were used to update the spacecraft attitude 
and gyro bias drift data, and to determine the orbital ephemeris of 
the spacecraft. The accuracy of the latter was found to be comparable 
with the accuracy of the ephemeris estimates ordinarily obtained from 
the data supplied by ground tracking stations. Thus it was concluded 
that the use of landmark and star' measurement data in this manner could 
provide two potential benefits for NASA's EOS program. Namely, an 
attractive alternative to ground tracking for the purpose of orbit de- 
termination and a significant improvement in the geometric accuracy cr 
fidelity of the telemetered earth imagery data. 

It was noted however, that the potential exploitation of ground 
truth data, for the purposes envisioned by the LUS study, imposed cer- 
tain requirements on the types of landmarks that would be most useful. 
Namely, those types that could be accurately and reliabily detected in 
the imagery data stream by means of an automated process involving a 
minimal amoiant of computer time and hardware. This latter requirement 
was due in part to the belief that several additional benefits could 
possibly be realized if the automated landmark detection scheme could 
be economically implemented in the EOS spacecraft itself rather than in 
a ground based data processing facility. These additional potential 
benefits relate to improvements in the precision of the onboard attitude 
control scheme and the fact that the location (or time of occurrence) of 
each landmark could then be conveniently recorded and subsequently iden- 
tified (in the data processing) by means of a specific word in the 
down link data stream. The possibility of realizing these additional 
benefits served therefore to augment the desire for an extremely rapid 
and efficient automated landmark detection scheme. 

Section 2.1 describes the historical reasons for believing that 
these stringent landmark detection requirements could be satisfied by 
using man made or artificial landmarks consisting of and defined by any 
sufficiently powerful source of radiation in the visible and/or near IR 
portions of the electromagnetic spectrum. In brief, it was theorized 
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that the electronic signature of such radiation sources (hereafter re- 
ferred to as beacon landmarks) should consist of a narrow (essentially 
single pixel) pulse, superimposed on the relatively low frequency 
(quasi d.c.) background signal generated by an EOS multispectral scanner. 
The' theorized pulse signature of beacon landmarks therefore provided a 
reasonably sound theoretical basis for envisioning a standardized auto- 
mated landmark detection scheme, based on the well known techniques of 
pulse detection, that could provide an intrinsically high precision 
with a minimal amount of computer time and hardware. It was realized, 
however, that the net value of a beacon landmark utilization program ■ 
would depend on the answers to a number of questions relating to the 
feasibility, practicality, and costs associated with the construction 
and performance of beacon landmark stations per se, and the implement- 
ation and operation of a national grid of approximately 100 such stations 
through the United States . The ALIS study was therefore undertaken with 
the general objective of generating answers to these types of questions, 
and the scope of the study was later expanded to include a series of 
beacon landmark experiments with Landsat 1 and Landsat 2 . The structure 
of the ALIS study and the conclusions reached are briefly summarized in 
the following section. 

1.2 . SUMMARY 


1.2.1 Beacon Landmark Trade Study 

, Section 2.2 describes the results of a trade study which examined 
the relative technical merits and costs of various candidate types of 
radiation sources for a beacon landmark station. The types of radiation 
sources considered by the study included visible and IR lasers, flat 
and curved mirrors, and high intensity searchlights covering the visible 
and near IR portions of the electromagnetic spectrum. The results of 
this study suggested that a Xenon searchlight would be the most prac- 
tical and cost effective type of radiation source for a beacon type 
landmark station. The practicality of using either visible or IR lasers 
as sources for a fixed beacon landmark was found to be severely limited 
by either their intrinsically low output powers or by their relatively 
low efficiences for converting electrical power into useful radiant 
power. In general, one pays a very high price per watt of radiant power 
for the highly monocromatic and collimated features of a laser beam, 
neither of which is actually required or even desired for a fixed beacon 
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landmark soxirce. The practicality of using a curved solar reflector 
as a source for a fixed beam artificial landmark was found to be severe- 
ly - limited by the uncommonly large surface that would be required to 
accommodate seasonal changes in the solar zenith angle. Smaller re- 
flector surfaces could be utilized, but only at the expense of a two 
degree of freedom sun tracking and pointing system. Thus, the -particular 
advantages that lasers and solar reflectors would otherwise offer, if 
the requirements called for a monochromatic and/or highly collimated 
beacon, turned out to be of little value for meeting the cost require- ' 
ments of a fixed beacon landmark. The size and pointing requirements 
for a Mirror Landmark Station or MLS are derived and discussed in detail 
in Appendix C. 

Subsequent investigations into the efficiency and maximum rad- 
iant power obtainable with conventional sources of broad band radiation 
revealed some rather surprising facts of particular value for the sel- 
ection of a cost effective landmark beacon. Throughout the past decade 
or so a significant technological advancement has been achieved in the 
area of concentrated high intensity light sources , primarily through 
the development of the xenon short-arc discharge projection lamp. Fea- 
turing an operating efficiency of nearly 50% , a very minimum of main- 
tenance over a long life (approximately 1000 hours) , and push button 
start up, xenon lamps can provide up to 10,000 watts of radiant power 
with a spectral distribution closely matching that of the sun. These 
desirable attributes of the xenon short-arc lamp account for the appear- 
ance, over the past decade, of many new types of optical equipment in- 
cluding new and compact types of military and commercial searchlights, 
signaling beacons, solar simulators, and wide screen projectors. Because 
of these factors it was determined (in Section 2.6) that a commercially 
available xenon searchlight, rated at 4 kilowatts (d.c. input power) , 
should adequately satisfy the intensity requirements for an EOS land- 
mark beacon, in the most practical and cost effective manner. The 
approximate cost of these xenon searchlights, including a power supply 
for operation from a 220 volt, 60 cycle single phase line, and a power 
relay or line contactor for remote control switching operations, is 
currently estimated at $7,000.00 (in lots of 100). 

1.2.2 Theoretical Signature Characteristics of Xenon 

Searchlight Landmarks 

Section 2.3 describes the theoretical model developed for pre- 
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dieting the pulse amplitude of a Xenon searchlight signature in the 
various output channels of a multispectral scanner (MSS) like the MSS 
in the current Landsat spacecraft and/or the MSS planned for the EOS 
spacecraft- For the purposes of this summary it is noted that the model 
transforms the actual (measured) bolometric intensity of a Xenon search- 
light into the theoretical amplitude of the pulse response expected in 
each MSS channel expressed as a percent of the full scale range for 
each, channel . • The pulse amplitude model is therefore formulated entire—, 
ly in terms of empirically known wavelength dependent parameters and in- 
tensity calibration data for the MSS and the Xenon searchlight under 
consideration . 

A significant prediction generated by this model is that the 
pulse amplitude in the near IR channel should be nearly twice as large 
as the pulse amplitude in the first (visible) channel. This contrasts 
significantly with the spectral response of pulses that may be readily 
generated by a variety of man made and natural specular reflectors of 
solar. radiation. In the latter case the near IR pulse amplitude would 
nominally be expected to be less than that observed in the first visible 
channel. This distinguishing characteristic of a Xenon searchlight land- 
mark is due to the existence of a pronounced XENON FLARE in the near IR 
portion of a high intensity Xenon discharge. Consequently, it was 
anticipated that this feature of a Xenon searchlight landmark could be 
exploited to improve the selectivity, and to therefore reduce the min- 
imum pulse amplitude threshold requirement for an automated detection 
scheme . 

1,2.3 Empirical Signature Characteristics of Xenon 

searchlight Landmarks 

The accuracy of the theoretical Xenon searchlight signature 
(pulse amplitude) model was determined by means of a series of Xenon 
searchlight experiments that were carried out with Landsat 1 and Landsat 
2 during the period from 3 November 1974 to 26 June 1975. The experi- 
mental program was conducted with the generous (unfunded) support of the 
New Hampshire Army National Guard (Searchlight Battery A, of the 197th 
Field Artillery Group) , which furnished the portable (tactical) Xenon 

searchlights and the necessary vehicles and volunteer personnel required 

. ( 2 ) 

for their transportation and operation. 
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The experimental searchlight landmark program resulted in the_ 
acquisition of 17 searchlight landmark data points in the Landsat MSS 
imagery data, providing 68 spectrally distinct searchlight signatures 
(in the four MSS bands) from 11 spatially resolved and non overlapping 
scan lines- This statistically significant data base was generated 
from eight separate experiments that were conducted at six different 
geographical locations (cities) encompassing the states of New Hampshire,. 
Massachusetts, New York, and Utah. The multiplicity of searchlight 
landmark data points resulted from the use of two or three spatially re- 
solved lights in a few of the experiments and a total of nine spatially 

(3^ 

resolved searchlights in the last experiment. Details concerning 

these experiments and the results obtained are presented in Section 2.4. 

1.2.4 Signature Modeling Accuracy 

The accuracy of the signature (pulse amplitude) model is de- 
termined in Section 2-5 by comparing the theoretically predicted pulse 
amplitudes in each wavelength channel of the Landsat MSS with the aver- 
aged values of the pulse amplitudes obtained from the digitized Landsat 
imagery data. On the basis of this analysis, the empirical pulse heights 
were found to be generally greater than the model predictions, with 
fractional .excesses of 7.7, 4. 8, -.03 and 10 percent (for the .55, .65, 

.75 and .95 micron channels) and corresponding standard ( la ) deviations 
of 26,22,21, and 15 percent. The decrease in the standard deviation' 
with wavelength was found to be consistent with the predicted and ob- 
served increase of pulse height with wavelength. Thus, to the accuracy 
indicated above, the xenon searchlight experiments with Landsat 1 and 
Landsat 2 served to establish the validity of the predictive model and 
the anticipated factor of -two enhancement of the band 7 pulse amplitude 
relative to the band 4 pulse amplitude. 

1.2.5 Selection of a Prototype Xenon Searchlight for an EOS SLS 

The above results were then utilized, in the manner described in' 
Section 2.6 to determine the size (or power rating) and cost of a Xenon 
searchlight that would be suitable for an EOS searchlight landmark sta- 
tion or SLS. In brief, the Landsat imagery data (obtained from the ex- 
perimental program) .was used to facilitate the design and testing of a 
digital Xenon Searchlight Landmark Detection (XSLD) algorithm. • The 
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XSLD algorithm, fianctions. as a filter, which incorporates a set of ab- 
solute and' relative (channel 7 to channel 6) pulse amplitude threshold, 
criteria, to identify the pixel numbers of all Xenon searchlight pulse 
signatures (if any) in the Landsat imagery data passed through the 
filter. The channel 7 pulse amplitude threshold was then varied- to 
determine the minimum' channel 7 pulse amplitude required to assure a 
high degree of discrimination against the variety of solar radiation 
pulse signatures found in the Landsat • imagery data. In this way it was 
found that a minimum channel 7 pulse threshold of seven percent (of the ' 
full scale) could be adopted as a guide for selecting a prototype. Xenon 
searchlight that would be suitable for an EOS SLS. 

Based on this result it was concluded that a commercially avail-, 
able Xenon searchlight (trade named XENOSPOT) would be suitable for an 
EOS SLS. In brief, it was found that the Xenospot has the theoretical 
capacity to generate a pulse amplitude of between 30 and 46 percent 
of the full scale range for the near IR channel of the MSS planned for 
the EOS spacecraft. The 30 percent value was derived from the conser- 
vative theoretical pulse amplitude model and actual measurements of the 
minimum calibration radiation intensity over a 1.5 by 6.0 degree rect- 
angular cone within the center portion of a slightly' defocused beam. 

These results were obtained however with a non optimal beam 
spread lens, and the manufacturer is confident that an optimally de- 
signed lens could boost the minimum intensity in such a cone by as 
much as 54 percent. The 46 percent estimate for the minimum pulse am- 
plitude over such a cone is therefore regarded as a realistic design 
goal estimate for the near IR pulse amplitude of an EOS SLS employing a 
Xenospot searchlight. The delivery price for a Xenospot, including its 
power supply and an optimally designed lens, would be approximately 7K 
in lots of 100 or more. 

The 1.5 by 6.0 degree cone size specification was adopted, in 
the above considerations, in order to provide a generous allowance for 
pointing offsets that can result from the anticipated cyclic deviations 
in the attitude and the orbital ephemeris of the EOS spacecraft. 
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Section 3 of this report attempts, to provide a relatively clear 
picture of the overall magnitude' of a searchlight landmark utilization 
program in terms of the major tasks that would have to be considered in 
connection with the implementation and operation of a national grid of 
searchlight landmark stations (SLSs) . The theory of operation and the 
basic types of components that might be incorporated in one possible 
configuration of an SLS are briefly described in order to provide a 
preliminary conceptual basis for the generation of an RFQ for (1) a 
detailed engineering design, (2) a. realistic cost estimate for the con- 
struction and testing of a prototype SLS, and (3) a realistic cost esti- 
mate for the fabrication, installation, and check out of 100 SLSs through- 
out the U.S.. Additional tasks relating to SLS site selection, land 
leasing, and surveying are also discussed. 
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SECTION 2 

SELECTION AND TESTING OF AN OPTIMAL 
. ARTIFICIAL LANDMARK FOR AN EOS SPACECRAFT 

2.1 MSS SIGNATURE CHARACTERISTICS OF BEACON TYPE ARTIFICIAL LANDMARKS 

2.1.1 Introduction 

This section reviews the general operation of a multispectral 
scanner (MSS) and the historically anticipated pulse signature character- 
istics of beacon type artificial landmarks which suggested their ideal 
suitability for the kind of landmark utilization program described in 
Section 1.1.1. An understanding of the geometrical relationships and 
imaging characteristics described here will be assumed in the later 
sections dealing with the theoretical pulse amplitude model, the search- 
light experiment results, and the discussions in Appendix A concerning 
the statistical probability of signature degradation due to detector 
sharing and data sampling. The imaging of extended sources (earth 
scenes) and the operation of a Landsat MSS responsible for the ordinary 
MSS scan signal is described first and followed by a description of the 
imaging of a beacon landmark., which results in the superposition of a 
pulse signature on the ordinary MSS scan signal. 

2.1.2 Extended Source Imaging Characteristics 

The various ray paths associated with the imaging of extended 
and point sources in a multispectral scanner (MSS) are ‘depicted by the 
simplified (heuristic) drawings in Figure 2-1. 

The Instantaneous Field of View or IFOV of an ,MSS is defined as 

the divergence angle of the individual detectors (or light pipes) 

located in the image plane of the MSS. The IFOV is therefore denoted 

by the subscripted angle parameter 0j^ (Theta detector) . The solid angle 

(or solid field of view) defined by each detector is then denoted by 

the corresponding solid angle parameter At each instant of 

• D 

txme a given MSS detector views a square element of the earth's sur- 

2 

face denoted by the area parameter Aj^q = R where R is the instant- 

aneous altitude of the spacecraft. A given detector is therefore 
illuminated by only that portion of the scattered solar radiation from 
Adg which is propagating and converging toward the MSS scan mirror with- 
in the directions defined by the viewing cone (as indicated by the 
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dashed lines (E) ) . The particular surface area (ADq) being viewed by 
a given detector in the iraage plane, at a given instant in time, is 
therefore determined by the instantaneous orbital position and attitude 
of the MSS, and by the instantaneous orientation of the scanning mirror 
{M) relative to the optic axis (OA) . For the Landsat MSS, R (mean) = 
.916 X 10® meters and = 86 x 10“® radians. Thus, = 7.4 x 10"^ 
steradians, and the instantaneous ground resolution width W = R0 q is 
equal to about 79 meters. 

The recorded portion of the electronic output from the MSS de- 
tectors is synchronized with that half of the scan mirror oscillation 
corresponding to an eastward sweep of the earth's surface {NW to SE in 
the crosstrack direction) , and the size of the strip thus scanned by 
each detector (in the Landsat MSS) is approximately 79 meters by 185 km 
(100 nra) . The varying illumination experienced by each MSS detector 
during a scan therefore generates an analog scan signal which is then 
amplified and digitized and telemetered to earth, where it is processed 
and used to generate a photgraphic scan line in the resultant imagery 
products. A representative analog scan signal from an MSS detector is 

illustrated in Figure 2-2 (following a calibration pulse) , and the com- 

4 

plete detector array for the Landsat MSS is illustrated in Figure 2-3. 

In reference to Figure 2-3, it is noted that the Landsat MSS 
employs six contiguous detectors in the intrack direction, which are 
split up or divided spectrally into four wavelength channels covering 
the following spectral ranges (in microns) : 

= 0.55 (+ .05) (2-1) 

X 2 = 0.65 (+ .05) 

X^ = 0.75 (+ .05) 

X^ = 0.95 (+ .15) 

Each oscillation of the scan mirror therefore results in the generation 
of six contiguous IFOV scan signals in each of the four wavelength 
channels, or a total of 24 spatially and spectrally independent scan 
signals per mirror scan. The motion of the scan mirror therefore 
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causes the total radiometfic image from a 0.468 by 185 ‘km (wide) strip 
to translate across the entire MSS image plane in the direction shown 
in Figure 2-1 (b) and 2-3. The scene illumination recorded in the 
successive wavelength channels (4 to 1) therefore corresponds to the 
illumination conditions at the scene for successively greater points in 
time. 


2.1.3 Beacon Landmark Imaging Characteristics 

The radiation arriving at the aperture of an MSS from a man- 
made beacon or solar reflector (as illustrated in Figure 2-1) can be 
described as plane parallel radiation from a distant point source. The 
ray paths of such radiation are indicated by the solid lines (P) . In a 
perfectly constructed optical instrument, the image of the source (point) 
would- be transformed into a series of concentric rings about a central 
spot (called an Airy disc) containing about 85 percent of the radiant 
power passing through the aperture lens,. The theoretical diameter of 
this central image spot is given by the well known formula 



2.44 X f 


"A 


( 2 - 2 ) 


where X is the wavelength of the radiation being considered; f is the 
focal length of the imaging system; and is the diameter of the aper- 
ture lens. In most practical optical instrinnents however, the image of 
a point source (or plane parallel radiation) is diffused throughout a 
larger spot (sometimes referred to as a blur circle) due to the combined 
effects of diffraction, figure imperfections, and misalignment of the 
optical elements. The diameter (D„) of the blur circle for the Landsat 

^ I 

MSS is approximately three times the diffraction value (D„) expressed 

5 ^ 

by Equation 2-2 . The ratio of the blur circle diamter (Dg) to the 

detector width (L = f 6^) for the Landsat MSS can be estimated from the 

formula 


a = °S 


7.62X 

°A 


(2-3) 
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Using the values of and 6^^ for the Landsat MSS (..23 meters and 86 
micro radians) gives the following values of a for each of the MSS 
channels ; 


“l. ^ 

0.21 

(2-4) 

“2 

0.25 


“3 = 

0.29 


“4 = 

0.37 



2.1. 3.1 Imaging Constraints 


The spot image of a landmark beacon is depicted in Figure 2-1. 
(b) at an instant when the beacon landmark .is located directly in the 
center of the viewing cone (and the instantaneous s.urface area , 

and when the image spot is, accordingly, positioned at the center of 
the detector. The observation of a point source of radiation, in this 
sense, is accomplished by means of the small fraction of its total 
radiation output which is propagating in the direction of the space- 
craft at the instant that the source point (or landmark) comes within 
the IFOV of a given detector. The radiation in all other directions is 
totally wasted at this point since it in no way contributes to the 
radiant power impinging on the scan mirror and thus the image spot at 
the detector. A certain minimal amount of beam divergence and thus 
wasted power would necessarily be required however in order to allow 
for the unavoidable uncertainties in the angular position and/or the 
attitude of the spacecraft relative to the beacon landmark at the instant 
of sighting. 

Thus for all practical purposes , the point source radiation 
impinging on the MSS, as indicated by the solid rays labeled (P) in 
Figure 2-1 (b) , may be viewed as being contained within an extremely 
narrow and directionally fixed beam, with a full cone angle (Dj^/R) on 
the order of 1.5 x 10'^ degrees. These remarks are made only to clari- 
fy a potential source of misunderstanding. In particular, it is noted 
that a beacon landmark cannot be repeatedly imaged or observed, on 
successive scans, unless the point on the earth's surface where it is 
located is repeatedly imaged on successive scans , and then only if the 
spacecraft is still within the macroscopic angular limits of the beacon. 
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2. 1.3. 2 Pulse Amplitude Considerations 


A pulse like signature for a beacon landmark would logically 
be expected to occur in the output of each MSS channel as the scanning 
mirror translates the corresponding image spot over the surfaces of the 
respective viewing detectors in the manner indicated in Figures- 2-2 and 
2-3. The amplitude of the pulse (in the ith wavelength channel) should 
then be directly proportional to the maximum amount of beacon power 

which is thereby superposed on the respective viewing detector'. 
This power can be expressed quite simply in terms of the aperture area 
= n D^/4) and the beacon radiation intensity I(l^) at the aperture. 
The value of I (Aj^) represents the average number of watts per square 
meter and per unit wavelength interval within the spectral bandwidth of 
the ith channel. 

Thus, 






BW. 

X 


(2-5) 


where k^ is the MSS optical transmission factor for the center wave- 
length A^, and BW^ is the corresponding bandwidth of either the radiation 
per se or the MSS channel filter, whichever is less. 


It is helpful to compare this expression with a similar express- 
ion for the background signal or extended source power which 

illuninates the full surface of the detector. 


P^oC^i^t) = k.RQ(A.,t)£Jj,A^ BW. ■ (2-6) 

Here, Rg(A^,t) is the average instantaneous radiance over the instantan- 
eous viewing' surface but as seen at the MSS aperture (that is, 

after being attenuated by the atmosphere). The value of Rg(A^,t) there- 
fore represents the intensity per lonit solid angle and per unit wave- 
length interval (watts/ (m -y -str) ) within the spectral bandwidth of 
the ith channel. The product of Rg(A^.,t) and ,may therefore be re- 
garded as an "effective" average intensity I^(A^,t) over the aperture 
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which is analogous to 1(1^^) except for the fact that it is not plane 
parallel and it serves therefore to illianinate the full surface of 
the detector element. 


The calibrated magnitude of R^(A.) associated with the maximum 
(full range or full scale) output for a given MSS channel is a known 
parameter of . an MSS which is here denoted by • Substituting 

Rgjjdi) for Rg(X^,t) in Equation 2-6 gives' the related expression 
for the detector illumination power ] associated with a full 

scale response in the ith channel of an MSS. 




k . 
1 






BW. 


(2.-7) 


The ratios (100 (100 can thus be conveniently 

formed and used to describe the ordinary MSS scan signal level SS(X.,t) 
and the superposed beacon landmark pulse amplitude BPAd^) as a fract- 
ional percent of the full scale outputs for the detectors in the ith 
wavelength channel. 


Thus, 


SS(X^t) = 100 




( 2 - 8 ) 


and 


BPA(Xj^) 


100 l(X^) 


(2-9) 


The pulse amplitude of a beacon landmark signature can therefore be ex- 
pected to be directly proportional . to the radiation intensity of the 
beacon [I(X^)1 and inversely proportional to either the solid IFOV 
(fijj) of the MSS (or the square of its linear IFOV Oj^) . The radiation 
intensity I*(X^) required to satisfy a given minimum pulse amplitude 
or design specification BPA*(X^) is therefore directly proportional to ■ 
fijj and given by the expression: 


I* up = Rgj^ U^) BPA*(X^)/100 (2-10) 

= .5 ftp 

Equations 2-9 and 2-10 are utilized and/or referred to extensively in 
the remaining subsections (2.2 through 2.6). 
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2 . 1 . 3 . 3 Pulse Width Considerations 


The width of a beacon landmark pulse affects the accuracy with 
which the geographical position of the landmark can be established in 
the MSS data stream, and it provides a useful parameter for the devel- 
opment of an appropriate detection scheme. The theoretically expected • 
pulse width (at half amplitude) is discussed accordingly in terms of ■ 
both an analog or continuous time base (seconds) and a quantized time 
base associated with a digital picture element (pixel) . 

The theoretical width of a beacon landmark pulse (at half amp- 
litude) should be very nearly equal to the so called Dwell Time for a 
detector (DT) , which is the time required for a point in the image to 
move across the surface of a detector, or the time for the IPOV to 
sweep across a given point on the earth's surface. Identifying -Tg as 
the time interval associated with the active (forward) portion of each 
scan, and ,9^ as the total angle scanned relative to the spacecraft, 
permits the dwell time (DT) to be expressed by the simple formula: 

DT = (8jj/6g) Tg. (2-11) 

From -Figure 2-4 it is clear that DT represents the time interval during 
which one half or more of the image spot, and thus the image power 
(X^) , remains on a given detector. The detector dwell time should 
therefore provide a reasonably good estimate for the half amplitude 
analog pulse width provided only that the detector-amplifier response 
time constant is significantly less than DT. For the Landsat MSS, 

= 0.2 radians (11.56 degrees), = 86 x lO"*® radians, and 
Tg = .0315 seconds. The corresponding beacon landmark pulse width, 
with respect to the analog output of the Landsat MSS,- would therefore 
be expected to be about 13.5 microseconds (at half amplitude). The 
position uncertainties associated with the analog and digital signatures 
of beacon landmarks are discussed in Section 2.1.4. 
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The digital signature of a beacon landmark, as it would appear 
in the telemetered data stream, can be predicted from the following 
considerations relating to the operation of the Landsat MSS- The 24 
detector array in the Landsat MSS, as shown in Figure 2-3, is sampled 
in the sequence indicated by the numbers to the right of each detector 
during a total sampling period denoted by the parameter Tgj^ which is 
approximately 9.95 microseconds. The relatively instantaneous sampling 
duration for each detector is therefore about 0.41 microseconds and 
each and every detector is therefore repeatedly sampled once every Tg^^ 
seconds. The spacing of the channel detectors in the scan direction is 
chosen accordingly to insure that a given scene moving across the image 
plane will be sampled in the same, position relative to the detectors in 
the successive wavelength channels. 

The time interval (Tgj^) between successive samples of a given de- 
tector determines the] distance (AX) that a given point in the imagery 
moves within or relative to' a detector detween successive samples . This 
distance can be expressed as Ax • = gL, where g = Tgj^/DT defines the 
fraction of the dwell time represented by sampling period. Thus, in the 
case of the Landsat MSS, DT = 13.5 ys and g = 0.74, so that AX = 0.74L. 
(The imagery in this case is said to be over sampled, with the IFOV for 
each sample or pixel element containing about 26 percent of the IFOV 
from the previous sample or pixel element) . 


For the purpose of analyzing the digital pulse signature of a. 
beacon landmark, it is convenient to parameterize the fraction' of the 
image spot falling on a given detector, during the jth sample, by the 
linear fraction FDj = Xj/Dg of its diameter (Dg) which overlaps a given 
detector relative to either its left or right edge. It is further con- 
venient to use the subscript value j = 1 to denote the first non zero • 
value PDjf and the parameter a=Dg/L <1.0, as before, to denote the 
relative spot diameter (in units of the detector width) . V7ith these de- 
finitions the sum of the first two image fractions takes the form: 


FD^ + FD2 


1 + 


(1-g) 

a 


( 2 - 12 ) 
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Thus, for the case of perfect sampling (neither over nor under sampling) 
g = 1, and there can only be (at most) two non zero samples of the spot 
image which are then consistent with the restriction FD 2 = 1 - The 

maximum value of FD^ is just g/a = 1/a, or unity, whichever is largest. 
Since the values of = aU^), from equations 2-4, range between 0.2 
and 0.4, it is clear that the perfect sampling of a spot image would 
generally produce two non zero samples of the spot consistent with the 
relation FD 2 = 1 - FD^ where the value of FD^ (and thus FD 2 also) is ran 
domly distributed over the interval from zero to unity. 

.However, for the over sampling case (where g = .74) we see that 
the sura of FD^ and FD 2 ' 

+ PD 2 = 1 + (2-13) 

will generally be greater than one, and for the smallest spot size , 
a^ = 0.21, would even be greater than two, implying that the full spot 
image and pulse amplitude could be sampled twi.ce. These considerations 
are illustrated in Figiure 2-4 for three spot sizes corresponding to 
a = 1, 1/3, and 1/5, and for a slightly smaller sample step size AX = 

5 

gL = 0.6L. The ordinate scale defines the fraction of the spot area 
associated with a given linear diameter fraction overlap (FD.=X./D ) , 

D J s 

where : 

fj (Aj) = I + i [Aj + Sin"^(Aj)] (2-14) 

and Aj = 2FDj-l 

Thus to the extent that the image spot intensity is uniform over 

the spot area A_ = ir D^/4, the spot area fraction f • (FD. ) also repre- 
o s 3 3 

sents the fraction of the maximum pulse amplitude BPA (described in 
Section 2.1. 3. 2) that would then be recorded in the jth sample or pixel. 
The beacon, landmark position uncertainties associated with the size of 
the IFOV and the digital sampling are discussed next. 


2-11 



FRACTION OP MAXIMUM PULSE AMPLITUDE (percent) 



. RELATIVE SAMPLE (PIXEL) NUMBER (]) 
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Shown. for 3 Fractional Spot Diameters (a == D^/L) and a Sampling 
Step Size. AX = 0.60L. ^ ^ 




2.1.4 Beacon Landmark Position Accuracy 


Measurement data relating to the orbit ephemerxs and the atti- 
tude of the spacecraft, as well as the operation of the MSS, can be 
used to compute the surface coordinates corresponding to the center of 
the instantaneous ground area associated with a given detector in 
the MSS at a given instant of time. With reference to one set of the 
channel detectors shown in Figure 2-3 (the set for example) we denote 


the relative locations of a given intrack detector by means of a row- 
subscript R. The computed longitude and latitude associated with the 
center of the IFOV for the detector at the time of its jth sample, 

are then denoted by the parameters .and ii' 4 j^(j) respectively, and 

the corresponding true values of the longitude and latitude are denoted 
by the parameters <f> 4 j^(j) and 'f* 4 |^(j)* "^he computation errors associated 
with the determination of the longitude and latitude for the .pixel P 4 j^Cj) 
may therefore be denoted symbolically as A<!» 4 j^ = 

= 'J'4R(j)-'('4R(j)- 


For the heuristic case of an exact polar orbit, the detector 
width (L) , and thus the instantaneous ground resolution distance (W) , 
serve to identify a specific longitude increment A 4 > 4 j^(j) in the scan or 
crosstrack direction, and a specific latitude increment ia the 

intract direction. Thus if the imagery of a specific (kth landmark)' 

m rp 

point on- the earth's surface, with coordinates and \li , is known to 
have- been somewhere within the detector's surface area at the 

instant of its jth sample, then the longitude and latitude of the center 
of the corresponding IFOV or the jth picture element can be measured 
(in the imagery processing) as = "fiR and to within a maxi- 

mum uncertainty of +.5 A^ 4 R(j) and ji*5^'!'4R( j) in the longitude and lati- 
tude. 


In general , however , it should be possible to make use of the 
finite spot size to determine the position of the image spot, relative 
to the observing detector,, to an accuracy of better than + 0.5L. For 
example, suppose the maximum pulse amplitude BPA is such that the over- 
lapping of the spot diameter on a given detection by a fraction FDj^ 
could be -readily observed as a partial pulse (f^^BPA) against the back- 
ground signal. Then one of two cases would then prevail with regard to 
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the position of the image spot in the intrack direction-. In one case, no 
portion of the pulse signature would be observed in the data from either 
of the adjacent intrack detectors and could con- 

clude that the center of the spot is therefore at least ( . 5-FDj^) Dg away 
from either of the intrack edges of The intrack or latitude uncer- 

tainty in this case would be reduced to +. 5L [l-2a { . 5-FDj^) ] . Thus if 
= 0.3, this uncertainty becomes +.5L(l-,4a) or+.4L and + . 3L for a 
equal to 0.5 and 1.0 respectively. 

■In the other intrack case, the image spot- would be overlapping 
an adjacent detector (say r+i^ least but by no more' than 

.5Dg, since that would just reverse the respective roles of the primary 
and adjacent detectors. In this case one would know that the 
center of the spot is within a distance of ( . 5-FD^)aL of the edge of 
the primary detector, or just 0.2aL (for FD^^ = 0.3). This suggests 
that the intrack or latitudinal position uncertainty can be minimized 
the most in the cases where the beacon landmark image is shared by two 
intrack detectors . 

The same types of arguments can be applied to the cross track 
or longitudinal uncertainties. Here, one would make use of the relative 
pulse amplitudes, making up the pulse signature profile or shape, to 
narrow down the position of the image spot relative to the detector at 
the time corresponding a particular sample, called the primary sample. 
The statistics associated the probability of intrack detector sharing 
are discussed in Appendix A. 

2.2 BEACON LANDMARK TRADE STUDY 

The objective of this task was to determine the relative merits, 
with respect to practicality and cost, of various candidate types of 
beacons for an EOS beacon landmark station. The specific types of 
beacons considered included solar radiation from specular reflectors,' 
lasers, and commercial searchlights. As a basis for comparison, the 
radiant power required from a given type of beacon was established in 
accordance with the following criteria. The radiation intensity pro- 
duced at an altitude of 700 km by a beacon should be sufficient to gen- 
erate a pulse amplitude signature of at least 50 percent of the full 
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scale range in one spectral channel of an EOS ’MSS having an instantan- 
eous field of view (IPOV) of 43 itiicroradians or less. The radiant 
power criterion then follows from the requirement that this minimum 
intensity criterion must be satisfied over a rectangular conic section 
or beam spread covering 1.5 by 6.0 degrees. These criteria served to 
establish the minimum total radiant power required from a beacon and 
thus, through other considerations relating to performance and effic- 
iency, the size and cost of a beacon that, would be suitable for an 
EOS landmark station. 

The 43 microradian IFOV anticipated for the EOS MSS is one 

half of that in the Landsat MSS. The radiation intensity required, in 

accordance with Eq. 2-10, is therefore one fourth of that which would 

be required for a Landsat beacon landmark. The fact that the intensity 

is specified at .a lower mean altitude, 700 km versus 916, means that the 

output radiation intensity (at a given reference distance) is decreased 

1 

again by a factor of (700/916) - 0.58. Thus the near field beacon in- 

tensity required for an EOS beacon landmark is only about 15 per cent 
of that which would be required for a Landsat beacon landmark. 

Because of the seemingly obvious economic and/or spectral 
attributes associated with solar reflectors and lasers , the feasibility 
and practicality of using these types of sources for a fixed beacon 
landmark were investigated first, with the following results. The 
practicality of using either visible or I-R lasers as sources for a 
fixed beacon landmark was found to be severely limited by ' either their 
intrinsically low output powers or by their relatively low efficiences 
for converting electrical power into useful radiant "power. In general, 
one pays a very high price per watt of radiant power for the highly 
monochromatic and collimated features of a laser beam, neither of which 
is actually required or even desired for a fixed beacon landmark source. 
The practicality of using a curved solar reflector as a source for a 
fixed beam artificial landmark was found to be severely limited by the 
uncommonly large surface that would be required to accommodate seasonal 
changes in the solar zenith angle. Smaller reflector surfaces could be 
utilized, but only at the expense of a two degree of freedom sun track- 
ing and pointing system. Thus, the particular advantages that lasers 
and solar reflectors would otherwise offer, if the requirements called 
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for a monochromatic and/or highly collimated beacon, turned out to be 
of little value for meeting the cost requirements of a fixed beacon 
landmark. The size and pointing requirements for a Mirror Landmark. 
Station or MLS are derived and discussed in detail in Appenxdix C. 

Subsequent investigations into the efficiency and maximum rad- 
iant power obtainable with conventional sources of broad band radiation 
revealed some rather surprising facts or particular value for the sel- 
ection of a cost effective landmark beacon. Throughout the past decade' 
or so a significant technological advancement has been achieved in the 
area of concentrated high intensity light sources , primarily through 
the development of the xenon short-arc discharge projection lamp. Fea- 
turing an operating efficiency of nearly 50%, a very minimum of main- 
tenance over a long life (approximately 1000 hours)., and push button 
start up, xenon lamps can provide up to 10,000 watts of radiant power 
with a spectral distribution closely matching that of the sun. These 
desirable attributes of the xenon short-arc lamp account for the appear- 
ance, over the past decade, of many new types of optical equipment in- 
cluding new and compact types of military and commercial searchlights, 
signaling beacons, solar simulators, and wide screen projectors. 

Because of these factors it was determined (in Section 2.6) that a 
commercially available xenon searchlight, rated at 4 kilowatts (d.c. 
input power) , should adequately satisfy the intensity requirements for 
an EOS landmark beacon, in the most practical and cost effective manner. 
The approximate cost of these xenon searchlights , including a power 
supply for operation from a 220 volt, 60 cycle single phase line, and a 
power relay or line contactor for remote control switching operations, 
is currently estimated at $7,000.00 (in lots of 100). 

2,3 THEORETICAL SIGNATURE CHARACTERISTICS OF XENON SEARCHLIGHT 

LANDMARKS 

This section describes the formulation of a theoretical model 
for predicting the maximum amplitude of the pulse that could be gener- 
ated in the ith wavelength channel of an orbiting multispectral scanner 
(MSS) due to the effective point source radiation from an xenon search- 
light beacon. The phrase "maximimn amplitude" refers to the pulse ampli- 
tude that would be expected for the case of a standard "clear" atmos- 


EM»R0X5UCIBILITY OF THE 
.OE-l’imAL PAGE IS POOR 



phere and the case where the image spot falls entirely within the in- 
track dimension (L) of one row of the intrack field detectors (as 
illustrated in Figure 2-3). Under these conditions, the beacon pulse 
amplitude 

BPA(X^) = (%FS) , (2-15) 

^D^GM ^^i^ 

as developed in Section 2. 1.3. 2, is expressed most simply as a fraction- 
al percent of the full scale intensity range for the ith wavelength 
channel- The only unknown quantity in the formulation at this point is 
the beacon radiation intensity I (Xj^) at the satellite. The solid field 
of view of the detectors (£5^^ the calibrated radiance ^^qj 5 (Xj^) 

corresponding to a full scale or maximum scan signal [SS (X^)=100%] , are 
accurately known parameters for a given MSS. In addition, it is noted 
that the radiance sensitivity implied by is determined primarily 

by the solar illumination conditions for a particular sun synchronous 
orbit. The values of for the Landsat MSS should therefore be 

essentially the same for all EOS MSSs which contain the same X^ channels 
and are designed for the same sun synchronous orbit and launch time 
(or window) . 

The basic objective of the present task is to simply express 
the beacon intensity at the satellite in terms of (a) its measured in- 
tensity Iq(Xj^) at some convenient and more or less arbitrary distance 
Rq from the searchlight, and (b) the average atmospheric optical thick- 
ness OT(Xji^) over the wavelength range or bandpass associated with the 
center wavelengths (Xj_) of the MSS channels. For this purpose it is 
convenient to express the MSS IFOV (Qq) in units of a reference IFOV 
0j3j^ = 25 microradians, by means of the normalized (dimensionless) IFOV 
parameter With this notation, the theoretical pulse 

amplitude model can be expressed quite simply as : 

BPAf = 1^(0) F^/ONjj)^ (%PS), (2-16) 

2 

where Iq( 0) [watts/meter ] is the measured bolometric (or spectrically 
integrated) intensity of the beacon at a distance from the search- 
light and at an angle 0 [degrees] from the central axis of the beacon. 
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In the case of a rectangular beam, 6 would be decomposed into two ortho- 
gonal components 0^^ and 02* The term contains all the other relevant 
constants and wavelength dependent parameters in the form; 


-OT. 
i 


[m^/watt] . 


BR 


GMi DR 


(2-17) 


Here, [1/microns] is the value of the normalized zenon spectral 

distribution (shown in Figure 2-5) ; exp (-OT^) is the atmospheric trans- 
mission due to a vertical optical thickness OT^ (averaged over BW^) ; 

BW^ [microns] is the MSS bandwidth for the ith channel? [meters] is 
the calibration distance associated with 1^(0); R[meters] is the alti- 
tude of the MSS above the earth's surface; BR_„. is the full scale 
radiance value [watts/m^-str] integrated (spectrally) over BW^; and 
OgRtstr] is the solid IFOV reference defined by 0^^ = (25xl0“® radians)^ 
The values of F^ and the corresponding parameters used in its comput- 
ation are listed in Table 2-1 for a reference calibration distance 
R = 62 meters, and a satellite altitude R = 916 km (which corresponds 
to the mean altitude of the current Landsat EOS spacecraft) . The opti- 
cal thicknesses (OT^) were obtained from standard atmospheric data 
reference tables^ 


A comparison of the values (hereafter referred to as the 
channel sensitivity constants) shows that the pulse height sensitivity 
increases monotonically with increasing wavelength, with a factor of 
nearly two in evidence between the .95 micron (or near IR channel sen- 
sitivity) and the .55 micron channel sensitivity. Inspection of Figure 
2-5 .shows that the relatively large sensitivity for the near IR channel 
is due to a pronounced "Xenon Flare" which falls completely within the 
.3 micron bandpass of the near IR channel of the MSS. The pulse ampli- 
tude sensitivity (PAS^) defined by: 

PAS^ E ^^^i = F^/ 0Ng (2-18) 

•y 

gives the pulse amplitude [%FS] per watt/meter of searchlight inten- 
sity at the reference calibration distance R^. For the Landsat MSS 
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NORMALIZED SPECTRAL DISTRIBUTIONS 
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TABLE 2-1 


VALUES OF PARAMETERS DETERMINING P. 

1 


CHANNEL (i) 

BW^[y] 



(D 

1 

o 

H- 

BRQjjiEW/M^- str] 

[m^/W] 

Eos 

Landsat 







1 

(4) 

.1 

.55 

1.05 

.685 

24.8 

2.118 

2 

(5) 

.1 

.65 

0.94 

.745 

20.0 

2.558 

3 

(6) 

.1 

.75 

0.88 

.785 

17.6 

2.867 

4 

(7) 

.3 

.95 

1.00 

.823 

46.0 

3,. 920 


R = 62 [m] R(Landsat) = 0.916 x 10® [m] = 25 x 10 ® [rad] 

O UK 

2i /0N^ ■< BPA. < 41 /0N^ 

o D 1 o D 





fijj = 86 micro radians , = 3.44 and PAS^ = 0.179, 0.216, 0.242, and 

0.331 respectively. The inverse (PAST^) = 5.587, 4.630, 4.132 and 
2 ^ 

3.021 [w/m*^/ {%PS) ] then gives the reference beam intensity associated 
with a theoretical pulse amplitude of one percent in each MSS channel. 

A pulse amplitude of 50 percent in the near IR channel, for example, 
would require a beam intensity (1^^) of about 151 watts/meter^ at a 
distance of 62 meters (203 feet) from the searchlight. Intensities of 
this magnitude are readily generated by military (tactical) xenon 
searchlights like the army's 30 inch (diameter) AN/TVS-3. Figure 2-6 
shows a family of curves representing the measured values of Iq( 9) 
from such a searchlight for various beam patterns parameterized by the 
number of turns (T) of a focus control. The data shown corresponds to 
a d.c. input power level of 15 kilowatts, which represents 75% of the 

maximum rated input power for an AN/TVS-3. The values of PAST^and 

-1 ^ 

PAS.y (given above) were used to construct the scales shown in Figure 

2-6 for the corresponding theoretical pulse amplitude response in 

channels 4 and 7 of the Landsat MSS. The experimentally determined 

accuracy of these scales, and the pulse amplitude model (defined by 

Equation 2-16) in general, is described in Sections 2.4 and 2.5. 

2.4 EXPERIMENTAL SIGNATURE CHARACTERISTICS OF XENON SEARCHLIGHT 
LANDMARKS 

In order to determine the accuracy of the theoretical pulse 
amplitude model represented by Equation 2-16, a series of AN/TVS-3 
searchlight experiments were planned and carried out with NASA's experi- 
mental EOS satellites called Landsat 1 and 2. The experiments were 
conducted with the generous (unfunded) support of the New Hampshire 
Army National Guard, Searchlight Battery A, of the 197th Field 
Artillery Group, which furnished all of the AN/TVS-3 searchlights and 

power supplies and the volvinteer support personnel required for their 

(21 

transportation and operation. 

2.4.1 Experimental Procedures 

The experimental procedures employed in the searchlight land- 
mark experiments and in the retrieval and location of the data are here 
briefly sianmarized as follows. Table 2-2 lists the satellite involved 
in each experiment, the date on which each experiment was conducted. 
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AN/TVS-3 SEARCHLIGHT INTENSITY AT 203 FEET 



THEORECTICAL PULSE HEIGHT {% FULL SCALE) 


£Z-\ 


TABLE 2-2 


landsat/searchlight landmark experiment data 


Landsat 

Date 

Test Site ' 

Searchlight Data 

CCT Data 


Number 

Mo/Dy/Yr 

Location 

No. 

Focus (T) 

Scene Number 

Tape 

Line/Pixel 

1 

11/03/74 

Manchester , NH 
W- Lon; 71.45“ 
N. Lat; 42.93° 

2 

3 (turns) 

1833-14491 

2 

1652/655 

2 

02/28/75 


1 

6 

2037-14494 

2 

1580/261 

1 

03/09/75 

W. Lon: 71.47“ 

1 

6 

1959-14435 

2 

1709/648 

2 

03/18/75 

N. Lat; 43.00° 

1 

6 

2055-14493 

2 

1587/254 

1 

05/03/75 

Westfield, MA 
W. Lon; 72.70° 
N. Lat: 42.16° 

3 

7 

5014-14472 

3 

489/543 

1 

05/07/75 

Batavia, NY 
W. Lon: 78.20° 

N. Lat: 43.02° 

3 

7 

5018-15094 

3 

1171/456 

2 

06/26/75 

Riverton , UT 
W. Lon; 111.95° 
N. Lat; 40.41° 

9 

9 

2155-17310 

3 

883/416 




the location of the test site for. each experiment, the number of search- 
lights employed in each experiment, and the searchlight beam spread 
or profile used in each experiment, as defined by the number of turns 
(T) of the. beam focus control (see Figure 2-6-)- i 

2 . 4 ,.l...l P.ointing 

In each experiment the searchlight (s) were aimed to intercept 
the predicted path of the Landsat orbit at right angles to the space- 
craft's instantaneous heading vector, which was very nearly equal to 
the point of closest approach to the test site. Predicted orbit posi- 
tion data furnished by NASA (GSFC) was used to determine the cross - 
track surface distance X between the test site and the subsatellite 
ground track. The elevation pointing angle (i|^) , measured relative to 
the local vertical, was then computed from the formula; 

sin(« = I (1+ [1+ (|)2 ( 1 + (2-19) 

where R is the estimated altitude of the spacecraft and is the 
radius of the earth at the latitude of the test site. The cross- 
track surface distance (X) was generally less than 30 km in all of the , 
experiments and the cross track orbit elevation pointing angle (4>) was 
generally between 1 and 2.5 degrees from the local vertical . The uncer- 
tainty in the true pointing angle was estimated to be within +0.1 
degrees . 


2. 4. 1.2 Firing 

Predicted orbit position data, pertaining to the time of the 
picture center coordinates for the Landsat scene covering each test 
site, was used to estimate the approximate fly-by time (t^) at which a 
test site would be scanned by the Landsat MSS. The values of t^ ranged 
from apprximately 10:40 a.m. in Manchester, N.H. to 11:20 a.m, in 
Riverton, Utah, and the searchlights were generally turned on from 5 to 
10 minutes prior to t^, and turned off about 5 minutes after t^. 
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2, 4. 1.3 Data Retrieval and Display 

Positive transparencies of the approximately 100 by 100 nauti- 
cal mile scene covering each test site were furnished- by NASA (GSFC) 
along with the corresponding system-corrected digital imagery data. 

The digital imagery- data was formatted on four computer compatible 
tapes (CCTs) with each tape containing the 25 n.m. wide imagery data 
from one fourth of the total scan. The complete data array for a scene 
picture, generated by one of the four wavelength channels, consisted of 
2,340 scan lines with approximately 3,240 intensity bits or pixels per 
scan line (810 pixels per tape). The pulse signature of 'each search- 
light was generally observable in the band 7 {near IR) transparency with 
the unaided eye; and careful measurements of its position, relative to 
the top and left edge of the scene, allowed its CCT line number and 
pixel number to be estimated to within about + 10 lines and +17 pixels. 
This process narrowed the search field considerably and permitted the 
actual locations of the signatures to be readily identified in the print- 
ed and plotted formats of the CCT data that were generated at CSDL. 

The identifying numbers for the Landsat scene and the CCT 
tape containing the searchlight signatures for each experiment are 
listed in Table 2-2 along liitith the approximate scan lines and pixel num- 
bers’ of the signatures. (The pixel numbers given are referenced -to the 
first pixel in the tape indicated) . Table 2-3 provides a sample -list- 
ing of the MSS data for the first experiment. The pixel intensities 
are expressed as a percent of the full scale output for the MSS channel , 
where the symbol $$ is used to denote saturation or 100 percent of full 
scale. Sample plots of the corresponding data for the target scan line 
and' the scan lines preceding and following it are shown in Figures 2-7 
to 2-9. The absence of any noticeable signal sharing by the adjacent 
intrack detectors (dr scan lines) is characteristic of all the data 
collected from the seven searchlight experiments. This data base which 
was expanded by the use of multiple searchlights included 17 spatially 
distinct signatures, or 68 spatially and spectrally distinct signatures, 
on a total' of 11 scan lines. 
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TABLE 2~3 


lANDSAT 1 MSS PIXEL INTENSITY DATA PROM TAPE 2 OF SCENE 1833-14491 
EXPRESSED IN PERCENT OF FULL SCALE OUTPUT (WITH $$ = 100) . 
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2.4.2 


Data Reduction and Analysis 


2. 4. 2.1 Definition of Experiigental Beacon Pulse Amplitude (EBPAf ) 

The experimental beacon pulse amplitude (EBPA^) observed in 
the ith channel of the MSS was consistently evaluated as follows'. Let 
P* represent the pixel number corresponding to the maximum signal level 
SL* of the searchlight signature in a given band of the MSS. The base 
line or background signal level BSL* was then defined as the average 
of SL(+P) and SL(-P)/ where SL{+P) was taken as the value of either 
SL(P*+1) or SL(P*+2), whichever was lower; and SL(-P) was taken as the 
value of either SL(P*-1) or SL(P*-2), whichever was lower. The empiri- 
cal pulse amplitude was therefore consistently defined by the equation: 

EBPA^ = SL* - BSL* = SL(P*) - . 5 [SL,(+P) +SL (-P) ] %FS (2-20) 

2. 4. 2. 2 Computation of Theoretical Beacon Pulse Amplitude (BPA^ ) 

The theoretical beacon pulse amplitude (BPA^) was computed 
using Equation 2-16 with = 86/25 = 3.44, the P^ values from Table 

2-1, and the intensity values = 1^(0, T) from the calibration inten- 
sity curves given in Figure 2-6. The calibration intensity for a 
given profile is directly proportional to the d.c. input power (PS) 
to the xenon discharge lamp. This input power# which is displayed on 
each AN/TVS-3 searchlight by means of current and voltage meters, was 
found to vary slightly from searchlight to searchlight and even from 
one firing, to the next of the same searchlight due to nominal deviations 
in the power supply control settings. The theoretical beacon pulse 
amplitude associated with a given searchlight was therefore adjusted in 
accordance with a power correction factor PCF = PS/15. The theoretical- 
ly expected pulse amplitudes were therefore obtained from the following 
formulas : 
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BPA 4 

== 0.179 

PCF 

ije,T) 

%FS 

( 2 - 21 ) 

BPA 5 

= 0.216 

PCF 

I ( 6 ,T) 
0 

%FS 

( 2 - 22 ) 

BPAg 

= 0.242 

PCF 

I (e,T) 
0 

%FS 

(2-23) 

BPA.^ 

= 0.331 

PCF 

i^(e,T) 

%FS 

( 2 - 24 ) 


The experiments and results obtained are briefly summarized in 
the following subsections. In those experiments involving more than 
one searchlight, it may help to recall that the pixel numbers correlate 
with a basically west to east scan. 

2.4.3 Experiment No. 1 (Landsat 1 ) 

The first experiment, of 11/3/74, was carried out at Grenier 
Air Field, Manchester, N.H., using two searchlights with PCF=1 and T=3. 
The searchlights were spaced 152 meters apart in the crosstrack dir- 
ection, and the signatures obtained are shown in Figure 2-8.* Due to 
a power, cable mishap, it was thought that the eastward light was turned 
on a few seconds too late to be observed. The results suggest that 
either this was true and the light output had only reached about 36 
percent of maximum, or the pointing angle was off by about one degree. • 
In either case no meaningful comparison could be made with theory." The 
westward" light, on the other hand, saturated all of the MSS channels. 
Consequently, the only conclusion that could be drawn from this data 
point was that the experimental pulse amplitudes were greater than 77 , 
81, 73, and 78 percent of full scale. This result is consistent with 
the predicted pulse amplitudes of 78, 94, 105, and 144 percent of full 
scale for 1^0=0, T«3) , but the exact amount of agreement "is indeter- 
minate. 


2.4.4 Experiments No. 2,3, and 4. (Landsats 2, 1,' and 2) 

Experiments 2, 3, and 4, of 2/28, 3/9, and 3/18/75, were ‘ 
carried out at the New Hampshire Army National Guard Armory in Manchester, 
New Hampshire, using one searchlight. The respective power correction 
factors (PCF) were 1.04, 1.04, and 1.03, and the beam spread for each 
light corresponded to the T =6 profile. The signatures obtained, with 
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the adjacent scan line data > are shown in' Figures 2-10 through 2-18. 
Experiments 2 and 4 were carried out with Landsat 2, and if one assumes 
that all the lights scored a direct hit, T=6) , then the 

theoretically expected MSS response is 46, 55, 62., and 85- -per cent of 
full scale for bands 4 to 7 respectively. A comparison of the theoreti 
cal and experimental responses is shown in Table 2-4 below, where- the 
symbols E and t are used to denote the experimental and theoretical 
pulse amplitudes (for convenience) and the symbol SR is used to denote 
the ratio E/T. 

TABLE 2-4 . 


Experimental Results and Theoretical Predictions • for Landsat 
Searchlight Experiments 2 , 3 , and 4 



(2/28/75) 


' (3/9/75) 




(3/18/75) 


EXE 

.No. 2 


. EXP. No. 3 


1 


EXP. No. 4. 


LANDSAT 2 


LANDSAT 

1' ■ 


■ - '• 

LANDSAT -2 - 



E 

T SR 

E T 

J 

SR 


E 

T SR 

Band 

4 

56 

46 1.22 

41 . , 46 

0_.89 


58 

45 , 1.29 

Band 

5 

73 

-55 1.33 

44 . 55 

0.80 


72 

55 1.31 

Band 

6 

74 

62 1.19 

, 45 62. 

0.73 


64 

61 1.05. 

Band 

7 

>90 

85 >1.06 

>83 85 

>0.,98 . 


>86 

84- ->1.02 



. - 



. 

4 


■ 


The Landsat 2 data was • observed' here to. be consistently higher than the 
Landsat 1 data by about 60 percent, on bands. 4 to 6 and by an. indeter- 
minate .amount on band 7 which saturated in all three experiments. . The 

BR (i) values for the Landsat 2 MSS are the same as for Landsat 1, and 

gm 

the same AN/TVS-3 searchlight was used in all three of these experi- 
ments. It can' only be concluded therefore that the searchlight inten- 
sity for experiment No. 3 was somehow erroneously reduced to about 66 
percent of its nominal value. . This could.be explained by an unbeliev- 
ably. large pointing offset, of about one degree, or by the assumption . . 
that the searchlight power cable was .inadyertently hooked up to the . 
lOkw transformer .receptacle instead of the 15kw receptacle.-' 
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2.4.5 Experiment No . 5 


Experiment No. 5, of 5/3/75, was carried out at Barnes Munici- 
pal Airport in Westfield, Massachusetts, utilizing three AN/TVS-3 
searchlights. The lights were set up at three corners (NE, SW, and SE) 
of a rectangle measuring 120 meters intrack by 157 meters cross track. 
The PCF values for the NE, SW, and SE searchlights were 1.00, 1.03, 
and 1.01, and the beam spread for each light was set at the T-7 profile 
The signatures obtained, and the adjacent scan line dataware shown in 
Figures 2-19 through 2-22. 

This experiment was carried out with Landsat 1, and at the 
time of the flyby the MSS was operating in the high gain mode for bands 

4 and 5, and with zero compression for bands 4, 5, and 6. The signa- 
ture for the NE light appears on scan line No. 488 where the band 4 and 

5 components are seen to be saturation limited because of the high gain 
mode employed. 

The signatures for the SW and SE lights appear on the next 
scan line (489) , and inspection of lines 487 to 490 reveals no visible 
evidence of intrack detector sharing. In this experiment, the intrack 
separation of the lights was equivalent to one and one half detector 
elements. Thus if DS is the effective diameter of the image spot and 
L is the width of the detector, this experiment shows that (1.5L + DS) 
is less than 2L, or that DS is less than .5L. This is consistent with 
the theoretical DS(X) .values (.2L, .25L, .29L, . 37L) computed in 
Section 2.1.3. If there had been an intervening line with no evidence 
of detector sharing, the same result (DS<.5L) would have been implied. 

A comparison of the experimental and theoretical pulse amplitudes is 
shown in Table 2-5 for bands 6 and 7 only. 
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.TABLE 2-5 


EXPERIMENTAL RESULTS AND THEORETICAL PREDICTIONS 
FOR LANDSAT/SE ARCHLIGHT EXPERIMENT NO. 5 

(May 3, 1975) 


NE LIGHT SE LIGHT SW LIGHT 




E 

T 

SR 

E 

T 

SR 

E 

T 

SR 

Band 

6 

22 

47 

0.47 

19. . 

48 0 

.40 

54 

47 

1.15 

Band 

7 

27 

64 

0.42 

29 

66 0 

.44 

>7 3 

64 

>1.14 



The 

signature of 

the SW light , 

as 

recorded by 

the 

Landsat 


MSS, is consistent with the two signatures recorded previously by the 
Landsat 2 MSS for experiments 2 and 4 , in that the experimental results 
are slightly greater than the theoretical prediction. This tends to 
confirm the belief that the scale factors or sensitivities are the 
same for the MSSs of Landsat 1 and Landsat 2. At the same time this 
suggests that the results of Experiment No. 3 were indeed abnormal, 
for reasons unknown. The experimental results for the NE and SE 
lights in this experiment are equally puzzling, however, because of the 
fact that they- are so identical and yet so far off from what was 
theoretically expected. Also, the actual band 6 and band 7 responses 
of these twO' lights are both very nearly equal to 38% of the actual 
band 6 and 7 responses for the SE light. It is also curious that the 
SR values for these eastward lights are very nearly the same as the SR 
values for the eastward light in the first experiment of 11/3/74. 

2.4.6 Experiment No. 6 (Landsat 1) 

Experiment No. 6, of 5/7/75, was carried out in the large park- 
ing lot of the Batavia Downs Race Track, in Batavia, N.Y., using the 
same experimental set up as that of Experiment No. 5. The signatures 
obtained, and the adjacent scan line data, are shown in Figures 2-23 
through 2-26. As in Experiment No. 5, there is no visible evidence of 
detector sharing, and the separate signatures occurred again on two 
consecutive lines, even though the searchlights were separated by the 
equivalent of 1.5 detector widths in the intrack direction. 
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In contrast to results of Experiment No. 5/ however, the signa 
tures of the two eastward lights came through loud and clear, and it 
was the signature of the SW light that was abnormally low, so low in 
fact that it- is not visible at all. Since power correction factors 
of 1.01 were reported for all three lights, it can only be concluded 
that an abnormally large pointing error was inadvertently introduced 
in the azimuth and/or elevation setting of the SW light. A comparison 
of the experimental results and theoretical predictions for the two 
eastward lights is shown in Table 2—6. 

TABLE 2-6 

EXPERIMENTAL RESULTS AND ''THEORETICAL PREDICTIONS 
FOR LANDS AT/SEARCHLIGHT EXPERIMENT NO. 6 

(May 7, 1975) 

NE LIGHT SE LIGHT 




E 

T 

SR 

E 

T 

SR 

Band 

4 

68 

35 

1.94 

51 

35 

1.46 

Band 

5 

>75 

42 

>1.79 

50 

42 

1.19 

Band 

6 

46 

47 

0.98 

45 

47 

0.96 

Band 

7 

>79 

64 

>1.23 

77 

64 

1.20 


These results show an abnormaly high response for the band 4 
and 5 signatures of the NE light and also the band 4 signature of the 
SE light. Inspection of scan lines 1170 to 1173 suggests that this 
can most likely be attributed to the unusually large amplitude and 
high frequency structure of the background noise. 

2.4.7 Experiment No. 7 (Lands at 2) 

Experiment No. 7, of 6/26/75, was conducted in the mountainous 
sagebrush terrain of the Army's Camp Williams Military Reservation in 
Riverton, Utah. This experiment utilized a total of 9 AN/TVS-3 search- 
lights and was carried out by the entire complement of the Airmy 
National Guard's Searchlight Battery A, in lieu of a night training 
mission previously scheduled for that date. The experimental 
arrangement employed is illustrated in Figure 2-27 an’d described below. 
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Figure 2-27 ARRANGEMENT OF SEARCHLIGHTS FOR LANDSAT 2 SEARCHLIGHT 

EXPERIMENT NO. 7 

In Figure 2-27 , the AN/TVS-3 searchlights are identified by the 
reference numbers 1-9. Lights 1-7 were spaced 157 meters apart in the 
scan direction and light No. 4 was aimed to intercept the-predicted 
orbit of Landsat 2. The beam spread of each light was adjusted to the 
T=9 profile, and the pointing angles for lights 5 to 7 and lights 3 to 
1 were increased and decreased symmetrically, relative to that of light 
No. 4, by 0.84, 1.4, and 2.25 degrees respectively. Lights 8 and 9 
were positioned along the intrack direction, spaced 134 meters {down- 
track) and 80 meters (uptrack) relative to light No. 3, and pointed in 
the same direction as light No. 3. The signatures obtained, and the 
adjacent scan line data, are shown in Figures 2-28 through 2-33. 

The first signature encountered, that for light No. 9, is 
shown in the data for scan line 882. The signatures for lights 1 
through 6 are clearly visible in the data for scan line 883, and the 
signature for light No. 8 is clearly visible in the data for scan line 
885. The 1^ (6=2.5, T=9) for lights 1 and 7 was very close to the end 
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points of the T=9 profile, where I^(6,T) drops off very rapidly- Thus, 
the absence of a signature for light No. 7 could easily be accounted 
for by a pointing offset (relative to the true position of Landsat 2) 
-of 0.3 degrees, inspection of scan lines 88-1 through 886 failed- to 
reveal any evidence of in track detector sharing. Thus, eleven scan 
lines in all, containing a total of 17 searchlight signatures-, were 
obtained wi-hhout any visible evidence of intrack detector sharing. A 
comparison of the experimental results and theoretical predictions for 
this experiment is shown in Table 2-7. 

TABLE 2-7 

EXPERIMENTAL RESULTS AND THEORETICAL PREDICTIONS 
FOR LANDS AT/SEARCHLIGHT EXPERIMENT NO. 7 
(June 26, 1975) 


Band 

Light No. 

1 

Light No. 

_2 

Light No. 

_3 


E 

T 

SR 

E 

T 

SR 

E 

T 

SR 

4 

.8 

6 

1.33 

11 

11 

1.00 

8 

15 

0.53 

5 

7 

8 

0.88 

15 

13 

1.15 

13 

18 

0.72 

6 

6 

9 

0.67 

16 

14 

1.14 

15 

21 

0.71 

7 

12 

12 

1.00 

18 

20 

■0.90 

31 

28 

1.11 


Light No. 

6 

Light No. 

5 

Light No. 

_4 


E 

T 

SR 

E 

T 

SR 

E 

T 

SR 

4 

12 

12 

1.00 

19 

16 

1.19 

19 

22 

0.86 

5 

16 

15 

1.07 

25 

20 

1.25 

21 

27 

0.78 

.6 

19 

16 

1.18 

29 

24 

1.21 

25 

31 

0.81 

7 

31 

.22 

1-41 

39 

30 

1.30 

43 

42 

1.02 


Light No. 

8 

Light No. 

_9 





E 

T 

SR 

B 

T 

SR 




4 

5 

14 

0.36 

6 

13 

0.46 




5 

11 

17 

0.65 

9 

15 

0.60 




6 

9 

19 

0.47 

12 

17 

0.71 




7 

15 

26 

0.54 

25 

24 

,1. 04 
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2.5 


COMPARISON OF THEORY AND EXPERIMENT 


2.5.1 Accuracy of Predictive Model 

The accuracy of the theoretical pulse amplitude model (Equation 
2-16) was estimated separately for each of the- MSS spectral bands (i) in 
terms of the average values obtained for the ratios 'of the experimental 
to theoretical pulse amplitudes (SRi) . The SRi values given in Tables 
2-4 through 2-7 are recompiled in Table 2-8 as a function of experiment 
number and the searchlight designation. The SRi values for the search- 
lights identified by an asterisk were excluded from the averaging process 
for one or more of the following reasons: (1) the exact value of the 
experimental pulse amplitude could not be determined because of satu- 
ration; as indicated by a ^ sign; (2) the CCT tape data was considered 
to be unreliable, as was the case for experiment number 5; (3) the ex- 
perimental pulse amplitude was either abnormally high or abnormally' low 
in one or more of the channels , based on an arbitrary 50 percent dis- 
crepancy criterion, as was the case for the east light of experiment 
number one, the NE light of experiment number six, and lights 7 to 9 of 
experiment number 7. There are a number of factors which could have 
accounted for such anomalously low or high results, but which could not 
be either conveniently controlled or identified after the fact. These 
factors include total or partial light failure, differences in the 
actual lQ(e,T) profiles between the many lights used, anomalous point- 
ing errors, an anomalously low or high contribution from the background 
noise (which was quite high for the hilly sagebrush terrain of experi- 
ment number 7) , and the possibility of MSS data collection and/or CCT 
conversion errors. 

The average SRi values of 1.08, 1.05, 0.97, and 1.10 shown in 
Table 2-8 were therefore based on the experimental signatures from 10 
xenon searchlights. The monotonic decrease in the corresponding stand- 
ard deviation (0.26, 0.22, 0.21, and 0.15) with increasing wavelength is 
consistent with the observed increase of pulse amplitude (and thus signal 
to noise) with increasing wavelength, as predicted by the theoretical 
pulse amplitude model. Part of the variation in the experimental data 
could also be attributed to undetermined variations in the actual verti- 
cal optical thicknesses of the atmosphere in one or more of the wave- 
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TABLE 2-8 


SUMMARY OF EXPERIMENTAL/THEORETICAL , 
RATIOS FOR SEARCHLIGHT PULSE AMPLITUDES 
(*Data Excluded From Estimates of Average Values) 


Exp . No . 

Light 

PCF 

SR 4 

SR 5 

SR 6 

SR 7 

1 

B* 

1.00 

- 

- 

- 

- 

1 

W * 

1.00 

> 0.99 

> 0.86 

> 0.74 

> 0.54 

2 

1 

1.04 

1.22 

1.33 

1.19 

1.06 

3 

1 

1.04 

0.89 

0.80 

0.73 

0.98 

4 

1 

1.03 

1.29 

1.31 

■ 1.05 

1.02 

5 

NE * 

1.00 

- 

- 

0.47 

0.42 

5 

SE * 

1.03 

- 

- 

0.40 

0.44 

5 

SW * 

1.01 

“ 

- 

1.15 

^ 1.14 

6 

NE * 

1.01 

1 . 94 ’ 

> 1.79 

0.98 

> 1.23 

6 

SE 

1.01 

1.46 

1.19 

0.96 

1.20 

7 

1 

0.94 

i .33 

0.88 

0.67 

1.00 

7 

2 

1.05 

1.00 

1.15 

1.14 

0.90 

7 

3 

1.00 

0.53 

0.72 

0.71 

1.11 

7 

4 

1.01 

0.86 

0.78 

0.81 

1.02 

7 

5 

1.07 

1.19 

1.25 

1.21 

1.30 

7 

6 

1.20 

1.00 

1.07 

1.18 

1.41 

7 

7 * 

1.00 

- 

- 

- 

- 

7 

8 * 

0.94 

0.36 

0.65 

0.47 

0.53 

7 

9 * 

0.84 

0.46 

0.60 

0.71 

1.04 

Averages 

: 

1.06 

1.08 

1.05 

0.97 

1.10 

STD DEV 

(+) 

0.06 

0.26 

0.22 

0.21 

• 0.15 
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length bands) from experiment to experiment. Any descrepancy between 
the actual and published values of any of these optical thicknesses / 
on the other hand, would be expected to show up as a fixed bias or 
•offset, which would be reflected in a consistently high or low value 
for the SRis and thus their averages . An additional source of variation 
in the data can also be attributed to a small but finite amount of 
2400 hz ripple that was present in the radiation from the AN/TVS-3 - 
searchlights . This should have been^ faithfully reproduced by the 
50,000 hz bandwidth of the MSS detector/amplifiers and resolved (or 
stopped) by the 0.4 .microsecond sampling interval (or shutter duration). 
The 2400 hz ripple results from the fact that the "d.c." discharge of 
the xenon lamp is sustained by an unfiltered full-wave rectified current 
from a 400 hz, three phase power supply. Measurements of the intensity 
ripple at 0-52, 0.7 and 0.9 microns revealed fractional variations, 
about the mean intensity of + .09, + .13 and + .08 respectively. Since 
the searchlight intensity variations associated with the ripple were 
statistically independent of the other sources of intensity variation, 
it can be argued that the contribution of such variations to the total 
standard deviation of the experimental results was less than 10 percent 
of the values obtained for channels 4 and 7 and less than 23 percent of 
the values obtained for channels 5 and 6. 

Thus, on the basis of the acceptable data from 10 xenon search- 
light signatures , it appears that the theoretical pulse amplitude model 
provides a satisfactory and slightly conservative basis for predicting 
the pulse amplitude of a xenon searchlight landmark in the near IR chan- 
nel of an orbiting MSS. Consequently it was decided to use the theor- 
etical pulse amplitude model, as it stands, for the purpose of select- 
ing a xenon searchlight that would be suitable for an EOS MSS. This 
selection is discussed next. 
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2.6 


SELECTION OF A LOW POWER XENON SEARCHLIGHT FOR EOS MISSIONS 


2.6.1 Calibration Intensity Considerations 

The present design 'specifications for the EOS MSS place an upper 
limit of 43 microradians on the IFOV of the first four MSS channels, 
with a design goal of 25 microradians. Apart from this, the first four 
wavelength channels of the EOS MSS should be optically equivalent to 
the four wavelength channels of the Landsat MSS, which have an IFOV of 
86 microradians. The bolometric calibration radiation intensity (I } 

. -I . . * ® 

required to assure a given minimum pulse amplitude (BPA^) in the near 
IR channel of the EOS MSS , can therefore be conservatively determined 
from the experimentally confirmed pulse amplitude model (Equation 2-16) 
as 


= BPA 4 * (eN^)^/P^, (2-25) 

* 

where is the minimum intensity over a rectangular 1.5 by 6.0 degree 

portion of the beam at a calibration distance (Rq) of 62 meters, and 
F^ is the channel constant for an MSS altitude of 700 km. Since the 
mean altitude of the EOS orbit is presently specified as 700 km, the 

I 

channel constant (F^) in Equation 2-25 is greater than that for the 
Landsat MSS by the factor (916/700)2=1.712. Inserting F .=1.712F. =6.712 , 
and (ONjj) =(43/25) =2.958 in Equation 2-25 gives the desired num erical 
relationship between the minimum calibration intensity 

1 ^* = .441 BPA* (2-26) 

and the minimum desired band 4 or near IR pulse amplitude. Comparison 
with Equation 2-24 shows that, for a given near IR pulse amplitude, the 
minimtim calibration intensity required for an EOS searchlight landmark 
is a factor of 6.8 times smaller than that which would be required for 
a Landsat searchlight landmark. The AN/TVS-3 searchlights that were 
used in the Landsat searchlight landmark experiments, at only 75 per- 
cent of their maximum power rating (15 kw vs. 20 kw (max)), were there- 
fore much too powerful, and hence costly, to be considered for an EOS 
searchlight landmark station (SLS) . 
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In order to determine the minimum calibration intensity actually 
required for an EOS searchlight landmark/ it was first necessary to 

'ft 

determine a realistic value for the minimum pulse' amplitude (BPA^) that 
could be reliably detected in the digital CCT data. For this purpose, 
a digital Xenon Searchlight Landmark Detection (XSLD) algorithm was de- 
veloped and tested using the CCT data obtained from the Lands at/Search- 
light experiments. 

2.6.2 Characteristics and Performance Results of an Automated • 

Xenon Searchlight Landmark Detection Scheme 

The XSLD algorithm performs basically as a xenon searchlight 
signature filter. The filter processes CCT data over a specified line- 
pixel data array, looking for narrow pulses, in all four bands, which 
exceed a set of specified amplitude thresholds and which comply with 
the known spectral characteristics of a xenon searchlight signature. 

The experimentally verified enhancement of the near IR pulse amplitude, 
due to the pronounced "Xenon Flare" {discussed in Section 2.3) was found 
to provide a highly effective basis for discriminating between low level 
xenon searchlight signatures and the variety of low level solar radia- 
tion pulse signatures that are present in the MSS imagery data. Experi- 
mentation with the CCT data showed that optimum selectivity and noise 
rejection was obtained by requiring that the channel 7 pulse amplitude 
be greater than 1.2 times the channel 6 pulse amplitude. 

The sensitivity of the filter was then determined by increasing 
the minimum pulse amplitude requirements (thresholds) until a total re- 
jection of erroneous pulses was achieved for at least 200,000 data 
points (pixels) associated with four 50,000 line-pixel data arrays taken 
from four separate and particularly noisy portions of the scenes from 
four different searchlight experiments. The pulse amplitude sensitivity 
of the filter, as defined by the resulting pulse amplitude threshold for 
the near IR channel, was thus found to be seven percent (of full scale) . 

Proper (positive) operation of the filter was verified during the 
selectivity tests by insuring that one of the four 50k data blocks used 
in the testing contained the eight xenon searchlight signatures generated 
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in experiment number 7. Positive filter operation was thus verified by 
observing that the filter automatically identified and flagged the 
correct scan line .numbers and pixel numbers corresponding to all 
eight of the xenon searchlight signatures contained in the data (see 
Figure 2-29, 2-30, and 2-32). The pixel number of a pulse, as such, 
was consistently defined by the maximum channel 7 pixel intensity, or 
in the unlikely event of a flat top pulse (consisting of two equal 
channel 7 pixel intensities) by the pixel number representing the first 
channel 7 pixel intensity maximum. 

2.6.3 Searchlight Selection 

Using BPA/ — 7 in Equation 2-26 gives a value of 3*1 (watts/ 
meter^) for the minimum calibration intensity (1^ ) of a xenon search- 
light that could be utilized in an EOS searchlight landmark station. 

From a practical point of view, however , allowance has to be made for 
an occasional 50 percent decrease in the recorded pulse amplitudes due 
to the phenomenon of detector sharing (discussed in Appendix A) with 
additional allowances for possible decreases in the beam intensity it- 
self, due to ageing and/or contamination of the searchlight optics, 
and/or possible increases in the optical thickness of the atmosphere. 

All in all, it would .not be unreasonable to consider a design calibration 
intensity specification of four times the absolute minimum (1^ ) , cor- 
responding to a nominal, theoretical, channel four pulse amplitude of 
28 percent of full scale. This requirement was found to be adequately 
satisfied by. a commercially available, 4 kw, xenon searchlight (trade- 
named Xenospot) . Calibration intensity measurements, made without the 
benefit of an optimally designed beam spread lens, yielded a minimum 

calibration intensity, over a 1.5 by 6.0 degree cone, of 13 watts/ 

2 * 

meter (4.2 ). Based on discussions with the manufacturer it appears 

that a xenospot calibration intensity design goal of 20 watts/meter^ 

(6.5 ) would not be at all unrealistic. Thus, from the foregoing 

theoretical considerations and experimental findings , it appears that 

the use of a xenospot xenon searchlight in an EOS SLS could be expected 

to generate a nominal landmark pulse of between 30 and 46 percent of 

the full scale range available to the near IR channel of the EOS MSS. 

This should be more than adequate to insure an extremely high degree 

( 12 ) 

of landmark detection reliability and station integrity. 
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SECTION 3 

PRELIMINARY CONSIDERATIONS RELATING 
TO THE DESIGN, IMPLEMENTATION, AND 
OPERATION OP BOS SEARCHLIGHT 
LANDMARK STATIONS (SLSs) 


3.1 INTRODUCTION 


The feasibility of generating readily detectable surface position 
marker pulses, in the multispectral scanner (MSS) imagery data from an 
earth observation satellite (EOS) , has been advanced theoretically and 
confirmed by means the Landsat/Search'light experiments reported in 
Section 21 The motivation for this work, as described in Section 1, 
was based on the findings of a Landmark Utilization Study which sug- 
gested that ground truth data of this nature can be exploited, in com- 
bination with star measurement data, to determine the orbital ephemeris 
of an EOS (thus easing the burden now placed on ground tracking stations) , 
and to improve the accuracy of both the EOS attitude control system and 
the geographic mapping process . The expected technical benefits and 
advantages of such a program need to be evaluated however in the light 
of the overall cost for the implementation and operation of a suitable 
network of EOS Searchlight Landmark Stations (SLSs) . This cost would 
of course include the cost for the initial design and constuction of a 
prototype SLS and its, qualification testing, as well as the subsequent 
generation of production engineering design drawings and finally the 
actual fabrication and installation of the total number of SLSs required 
for the landmark grid. Contractors are understandably reluctant to offer 
even as much .as a guess as to what they think the total cost of such 
vaguely defined tasks might be. This section of the ALIS report is there-- 
fore concerned not with the generation of realistic cost estimates, but 
with the generation of a more detailed description of an SLS and the 
specific tasks that would have to be undertaken in connection with the 
implementation of a useful SLS network for the EOS program. The approxi- 
mate costs of components and tasks which can be reliably estimated are 
noted, but it will be clear that they represent only a small part of the 
total implementation cost. It was originally hoped that this phase of the 
ALIS study would result in a rough, first cut, ball park estimate for 
the overall implementation cost, but it has become increasingly clear 
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that this can only be reasonably accomplished through the generation of 
(and the responses to) RFQs for the major tasks. The purpose of this 
phase of the ALIS study is therefore to clarify the specifics of those 
major tasks to which such RFQs would be addressed, and in doing so to 
provide a more detailed assessment of the total magnitude of an SLS 
implementation program. 

An assessment of the approximate distribution and the number of 
SLSs that would be required for a landmark utilization program is gen- 
erated in Section 3.2. Section 3.3 describes the presently envisioned 
theory of operation of a single SLS in terms of a preliminary concept- 
ual reference design, and Section 3.4 discusses the practical problems 
of SLS site selection and the acquisition of land and/or leases for 
such. Section 4 then summarizes the conclusions that have been reached 
in the ALIS study and the recommendations for future studies and RFQs 
that are needed to obtain realistic cost estimates for the overall im- 
plementation task. 

3.2 ASSESSMENT OP SLS GRID SIZE 

f 1 1 

The analytical results obtained in the Landmark Utilization Study 
have indicated that the locations of the landmarks could be restricted, 
for the most part, to the continental U.S., with a frequency distrib- 
ution sufficient to insure the observability of only one or two land- 
marks per orbital pass. A seemingly reasonable SLS distribution, is 
illustrated in Figure 3-1 in which three SLSs are distributed along 
most of the 42 EOS ground tracks in the U.S., to provide a nominal 
redundancy of between one and two sightings per orbital pass over the 
U.S., using a total of 114 SLSs. A redundancy of two sightings could 
very well be excessive, in many instances, and future LUS studies would 
have to be undertaken to determine the statistical variations of the 
cloud cover along the EOS ground tracks, and also, the minimum number 
of sightings per orbit that could actually be tolerated for the kinds 
of applications and accuracies desired. In view of these considerations, 
it might be reasonable to develop fabrication and installation oriented 
RFQs on the assumption that 100 SLSs would be required. 
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Figure 3-1 Approximate Distribution of 114 Searchlight 

Landmark Stations. Slashed Lines Indicate 
■ Directions of EOS Ground Tracks. 
























3.3 SLS CONCEPTUAL REFERENCE DESIGN 


This Section provides a qualitative description of the presently 
envisioned construction and operation of an SLS, consistent with the 
operational guidelines discussed in Section 3.3.1. The description, as 
such, provides only a preliminary SLS conceptual reference design, 
which could possibly serve as the basis for generating an RFP to obtain 
a detailed engineering design and realistic cost estimates for the con- 
struction of a prototype SLS and also for the fabrication, installation 
and check out of 100 production models. 

3.3.1 Operational Guidelines 

For the purposes of this study, it was assumed that the EOS satel- 
'lite' design will include provisions for the continuous transmission 
(over the u.S.A. at least) of an omni-directional carrier signal, which 
can be employed to activate each SLS one or two minutes before the 
satellite passes over the SLS. Activation, in this sense, means the 
initiation of an automated timing sequence for controlling the start- 
up, operation, and shut-down of the SLS. It was also assumed that the 
SLS design would include an optional provision for communicating its 
identify or station number and its operational status or integrity to 
an onboard data collection system (DCS) , like that currently employed 

in the Landsat spacecraft for collecting and disseminating data from a 

12 

variety of earth based data collection platforms (DCPs). In addition, 
it was assumed that the SLS design would provide an adequate degree of 
protection against potentially harsh environmental factors, including 
provisions for an automatic NO-GO or abort in the case that inclement 
weather conditions prevail at the time of the flyby. 

The basic operating principles and components associated with the 
presently envisioned SLS design are briefly described next with refer- 
ence to the function block diagram shown in Figure 3-2. 

3.3.2 SLS Power Requirements 

Primary station power is provided by a 220 volt single phase line 
capable of furnishing a peak power load of approximately seven kilowatts 
for a period of about six minutes once every 18 days when the search- 
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figure 3-2 Searchlight Landmark Station Block Diagram 




















light is turned on. An internal transformer (not explicitly shown) is 
used to provide a 500 watt source of 110 volt power for SLS activation 
and control components. With the exception of a 14 minute standby and 
operation period, which occurs only once every 18 days* the SLS power 
consumption should be quite small , being only that which would be re- 
quired to operate a 60 Hz station clock and to maintain a full charge 
on the clock-drive back-up battery. The power, if any, that would be 
consumed by the SLS heater is difficult to estimate at this time, since 
it will depend on both the actual location of the SLS and the operating . 
temperature ranges for its components . 

3.3.3 SLS Activation Window 


In accordance with the guidelines given in Section 3.3.1, the 
activation of each SLS in advance of an EOS flyby will be triggered by 
an omni-directional carrier signal emanating from the EOS . Thus , in 
order to limit the activation of SLSs during each orbit to only those 
which lie within the actual ground path of the EOS, and to also maxi- 
mize the overall efficiency and lifetime of each SLS, the SLS design 
should include a timing mechanism that prevents electrical power from 
being applied to the SLS RECEIVER (and all other SLS components) unless 
the local date and time at an SLS coincide with that for an EOS flyby 
to within zero days and approximately +5.0 minutes. 

A number of possible configurations can be envisoned fof~ such a 
timing mechanism, and there is a natural tendency to envision a con- 
figuration based on the use of modern digital hardware, encompassing 
a precision digital clock, two preset counters, a permanent memory 
storage buffer, switching relays and associated d.c. amplifiers, and 
a back-up battery sufficient to supply the power that would be required 
to keep all timing components functioning during an extended power 
failure. The specific problem at hand is considerably simplified how- 
ever 'because of the precise 18 day periodicity of successive flybys and ^ 
the fact that the variations in the local flyby times are controlled 
and kept within a small fraction of the +5.0 minute activation window? 

*These considerations would not 'be seriously affected by a daily 
activation for 18 EOSs. See discussion in 3.3.4. 
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Thus, while an activation window of this size represents a very small 
fraction of 18 days, it also represents a very generous margin of safe- 
ty for the nominal timing errors that can be associated with a simple 
and reliable mechanical gating mechanism, driven by a synchronous 60 
Hz electric clock motor. 

The basic .components and operation of such a timer are illustrated 

in Figure 3-2 . A primary 60 Hz station clock is used to drive a pair 

of mechanically coupled conincidence gates or switches. The 'ORBIT GATE' 

switch, being driven at one revolution per 18 days, closes for a period 

of about 10 hours, centered on the nominal flyby time {T^) , T^ refers 

o o 

to the ith satellite in a pattern consisting of 2, 3, 6, 9, or 18 EOSs 
(see discussion at the end of 3.3.4). When this happens, power to the 
station receiver is then controlled exclusively by the TIME GATE switch, 
which, being driven at the rate of one (or two) revolutions per day, 
closes accurately for a period of about 10 minutes, also centered on 
When this happens, line power is then applied to the station RECEIVER 
and the SLS is said to be in a STAND-BY condition. The receiver power 
on/off switches are in parallel with the Activation Gate switches shown 
in Figure 3-2. Integrity of the station clock during a local power 
failure is secured by a back-up, battery powered oscillator and a 60 Hz 
amplifier. 

3.3.4 SLS Activation Functions 


Detection of the EOS carrier transmission can be said to occur at 
that instant when the strength of the received signal exceeds a certain 
threshold value, causing closure of a relay in the receiver and a low 
voltage d.c. activation signal to appear at the receiver output (R) . It 
is envisioned that the receiver gain would be adjusted so that detection 
will occur, nominally, when the EOS is approximately two minutes or 
about 890 curvilinear miles away from an SLS. Following detection of 
the EOS carrier signal, the activation signal is applied to a MOSITURE 
TRANSDUCER, designed to block passage of the activation signal if there 
is a detectable amount of snow and/or rainfall present. Otherwise, the 
activation signal forces closure of a normally open WEATHER INTERRUPT 
SOLENOID SWITCH which feeds a low current 110 volt signal (through the 
now closed activation switches in the Orbit and Time gate components) 
to activate the six minute station control timer. The circuitry in 
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this TIMER should be such that once the TIMER is energized it provides 
full and independent control of the remaining SLS functions , causing 
line power to be applied to the control timer itself/ the searchlight 
power relay, the dome-hatch drive motor, and the SLS Data Collection 
Platform components . A Maximum Signal Discrimator circuit is envision- 
ed for restricting the transmission of station integrity data until the 
BOS is approximately overhead. The SLS response time, for searchlight 
warm-up and dome-hatch opening, is estimated at 0.5 minutes. Con- 
sequently, the combined timing errors associated with the SLS clock, 
the time gate switch, and the EOS flyby can be as large as + 9.5 min- 
utes and as small as - 4.5 minutes. Nominally, the time gate cam would 
be offset by two minutes so that the total timing error range is sym- 
metric (+7.0 minutes). 

Allowance for multiple EOSs, in a symmetric pattern consisting of 
2, 3, 6, 9, or 18 EOSs could easily be provided by incorporating 18 in- 
dividual but adjustable risers or cams in the "18 day" cam that could be 
easily set, by means of a simple twist-lock action, in either the "out" 
(activate) or "in" position. The frequency of the standby condition 
(or activiation window) could thus be manually increased as needed, or 
as anticipated, to allow for successive flybys every 9, 6, 3, 2 or 1 
day(s). In this regard it would be advisable to make the stand-by 
window equal to the launch window. 

3.3.5 SLS Housing 

The SLS housing envisioned and illustrated in Figures 3-3 to . 3-9 
is briefly discussed as follows. Figure 3-3 shows an artists conception 
of an SLS located in a field. In practice, one half or more of the SLSs 
would probably be located on top of commercial buildings . The estimated 
dimensions for the overall SLS and the individual subsections is .shown 
in Figure 3-4 . The alignment baseplate shown in Figure 3-5 is prealign- 
ed so that when the SLS is installed, with the *'d leg on the *'d raoimt- 
ing bolt, the 6 degree beam spread from the xenospot will be automatical- 
ly aligned in the cross-track direction for that SLS site. The vertical 
dimension of the searchlight power supply section shown in Figure 3-6 
was configured initially to accommodate the standard three phase power 
■supply which is required for continuous operation of the searchlight. 

The manufacturer has stated, however, that a smaller (20 inch high). 
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220 volt single phase power supply would be adequate for the SLS duty 
cycle. Figure 3-7 shows the approximate physical arrangement of'th,e 
electronic components shown in the block diagram of Figure 3-2, and 
Figure 3-8 shows the approximate layout of the. housing section contain- 
ing the 4 kw xenospot searchlight. Figure 3-9 then shows an approxi- 
mate layout of the dome-hatch section with the BOS receiving and data 
collection platform transmission antennas. During operation the dome- 
hatch is rotated 180 degrees to a full open position as indicated in 
Figure 3-8. 

3.4 SLS INSTALLATION AND OPERATION COST CONSIDERATIONS 

The largest expense associated with the installation and operation 
of a national SLS grid would, in all likelihood, be that associated 
with the purchase and/or leasing of the land parcels and/or building 
spaces on which the SLSs would be located. A realistic estimate for 
these fixed and yearly costs could be obtained through a study encom- 
passing the following two work statements: 

(1) Assiuning that 114 SLS sites will be required, and 
that the down range or intrack separation between SLS sites 
along any EOS ground track must be greater than 300 miles, and 
that the offset of any SLS site in the crosstrack direction 
must be within + 10 nautical miles of the nominal ground track, 
determine the optional (most cost effective) distribution of 
SLS sites with regard to a realistic assessment of the RELATIVE 
(as opposed to absolute) purchase costs and leasing options for 
the various candidate types of SLS real estate within the total 
SLS location zone defined by the intrack and crosstrack con- 
straints noted above. 

(2) From the most effective distribution of SLS sites de- 
termined in (1) which is based on RELATIVE real estate purchase 
and leasing costs , select a representative sample of 15 potential 
SLS sites and, through actual negotiations with the property 
owners, determine the actual purchase and leasing costs for the 
15 sites selected, and then weigh these results in an appropriate 
manner to obtain a realistic (+ 10 percent) estimate for the 
total purchase and leasing costs for the entire grid of 114 SLS 
sites . 



Additional factors that would have to be included in both of the 
above tasks would be (a) the cost of bringing in a 220 volt, single 
phase power line to the test sites, (b) the yearly electric power costs, 
(c) the accessability of the site locations to major roads (for the 
purposes of installation, repair, and surveying) , (d) the costs for 
access right-of-ways where necessary, and (e) the total cost for an 
accurate survey of the 114 SLS locations . 
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Figure 3- 3 SLS Installation 
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Figure 3-4 SLS Dimensions (1/12 Full Scale) 
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Figure 3-6 SLS - Section A (1/8 Full Scale) 
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Figure 3-8 SLS -Section C {1/8 Full Scale) 
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SECTION 4 


CONCLUSIONS AND RECOMMENDATIONS 


4.1 CONCLUSIONS 

The feasibility of generating readily detectable searchlight 
landmark position marker pulses in the multispectral scanner imagery 
data from an earth observation satellite has been advanced theoretically 
and confirmed by means of a series of experiments with Landsat 1 and 
Landsat 2. And the results were extrapolated to identify a commercially 
available low power xenon searchlight that would be suitable for an EOS 
Searchlight Landmark Station (SLS) , based on an EOS MSS IFOV of 43 
micro radians. 

A n\imber of the major tasks that would have to be undertaken in 
order to carry out the actual implementation of a suitable national 
grid of SLS were discussed in terms of a preliminary conceptual refer- 
ence design for an SLS, and an attempt was made to outline the types 
of future studies that would have to be carried out in order to obtain 
a realistic assessment of the approximate costs that would be involved 
in such a program. These studies are summarized next. 

4 . 2 RECOMMENDATIONS 


The technical benefits and advantages that the EOS program could 
expect to derive from the type of Landmark Utilization Program discussed 
in this report need to be carefully weighed against a realistic estimate 
for the total expenses that could be expected in connection with the 
implementation and operation of a national grid of SLS. It is perhaps 
regretable that such an estimate can not be reliably obtained without 
further studies, but the nature of the problem is such that realistic 
estimates simply cannot be generated when such a large gap exists 
between the general concepts that presently exist and the specifics 
that have to be identified and properly appraised at a point much closer 
to the potential suppliers of the hardware products, services, and real 
estate property that would actually be required. In this regard the 
following studies are recommended for future consideration. 
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(1) An extenstion of the Landmark Utilization Study 
to establish a realistic distribution for SLSs throughout 
the continental U.S. taking into consideration such factors 
as cloud cover statistics and the corresponding minimum 
SLS redundancy factor. 

(2) An RFQ for the preliminary design of a prototype 
SLS with realistic {+ 10 percent) cost estimates for its 
fabrication and qualification testing, and for the gener- 
ation of engineering production design drawings and the 
fabrication and installation of the quantity of SLSs 
determined from (1) above. 

(3) An RFQ for a study to determine a realistic esti- 
mate for the costs that would be associated with the 
purchasing and/or leasing of SLS real estate, electric 
power hook up, and surveying, as detailed, in Section 3.4. 

(4) An RFQ for the design and testing of a prototype 
near real time xenon searchlight pulse signature detector 
for use in an on board detection scheme. 


^m^RODUCIBILlTY OF THE 
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APPENDIX A 


STOCHASTIC DEGRADATION OF 
BEACON LANDMARK PULSES 
DUE TO DETECTOR SHARING 
AND DATA SAMPLING 

A.l DETECTOR SHARING STATISTICS 

The maximum amplitude of a given beacon landmark pulse should 
occasionally be degraded by as much as 50 per cent due to the fact that 
the corresponding power (or number of photons per second) will occasion- 
ally be divided- equally between the spectral detectors in two adjacent 
rows of the totai intrack-field/spectral detector array as illustrated 

i 

in Figure A-1. Referring to Figure A-1, -we note that v?hen a beacon 
landmark is sighted, the position of its spot image along the in track 
coordinate of the MSS detector array is determined entirely by the 
position which the landmark happens to occupy with respect to the total 
intrack field of view (9ip) at the instant it is sighted. Thus, because 
of the Stochastic variations in the orbital ephemeris and attitude 
of' the EOS, it is reasonable to assume that artificial landmark images 
will have an equal a priori probability of occurring anywhere along 
the intrack coordinate of the MSS detector array. 

The statistics of the pulse amplitude degradation is derived here 

for the realitic case in which the maximum spot diameter (D^) is less 

than the linear dimension (L) of each detector element. The analysis 

is facilitiated by assuming a uniform photon flux density over the 
2 

spot area A„ = ttD /4. With reference to Figure A-1, the parameter FD 
s s 

(expressed in units of D^) represents the linear distance by which the 
edge of an image spot protrudes into any given detector row. Since the 
occupation space of an in-track detector row is L, the probability that 
the edge of an image spot overlaps by a distance between FDD^ and D^ 
from either edge is numerically equal to the probability that the cen- 
ter of the spot (regarded as point) is located within an enclosed path 
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of width S = [L - Dg(2 FD-1) ] • Since every possible position of the 
spot along the intrack coordinate is equally probable, and the prob- 
ability of finding the center of the spot in some detector row is unity, 
it follows that the probability of finding the center of the spot in a 
strip of width. S is just equal to the ratio S/L. Thus, 

P (D^, PD) = S = [1 - D_ (2 FD-1) 3 (A-1) 

“ t " —B. 

Equation (A-1) shows that the probability of finding the entire spot 
in an in-track row of the detector array is simply 

P (D^, PD=1) = 1 - (A-2) 

while the probability P (D , FD = .5) of finding from 50 to 100 percent 

of the spot in a detector row is clearly always equal to unity or 100 

percent. FD) can therefore be used to compute the probability 

that the beacon landmark power falling on a detector row will exceed 

a certain fraction (RPH) of the total power in the spot image, defined 

as the Relative Pulse Height (RPH) . For the assumed case of a uniform 

flux density over the spot, RPH is simply the fraction of the total 

2 

spot area (A = ir D /4) which is contained in a spectral detector row 

5 S 

due to the protrusion of the edge of the spot by a distance D^FD. RPH 
is therefore equal to FD only for the special cases when FD is equal 
to either 1-0 or 0.5. With the assumption of uniform intensity then, 
the relationship between RPH and FD is given by the equation 

RPH(FD) =7+1 [A / 1-A^ + sin“^(A)] (A-3) 

where A = (2 FD-1) ■ (A-4) 

Thus, if FD' and RPH' (FD') are the actual sampled values of .tbe 
fractional overlap for the spot diameter and area (or power) respect- 
ively, then the probability P(D ,FD) given by Equation (A-1) is the 
probability that FD' is between FD and 1.0, and the probability that 
RPH' is between RPH (FD) and 1.0. The quantity 1-P then gives the 
probability that the observed amplitude of a beacon landmark pulse 
(BPA' ) was degraded by a factor greater than l-RPH(FD). 



Table A-1 lists six values of FD and RPH(FD) in the range of 0,5 
to 1.0 for the smallest and largest spot dismeters of the. Lands at MSS, 
as given by Equation 2-4 of Section 2. Specifically, ( . 55u)/L=0, 21 
and DS(.95vi)/L = 0. 37- 


TABLE A-1 

STATISTICAL DEGRADATION OF ARTIFICIAL LANDMARK PULSES 
DUE TO THE PHENOMENON OF DETECTOR SHARING 


Linear Overlap 
Fraction 

Probability that observed pulse 
height is between RPH and 100% 
of its nominal (undegraded) value 

Relative 

Pulse 

Height 

FD 

P(.55u) (%) 

P(.95y) (%) 

RPH(FD) (%) 

1.0 

• 79 

63 

100 

0.9 

83 

70 

95 

0.8 

87 

78 

86 

0.7 

92 

85 

75 

0.6 

96 

93 

63 

^ 0.5 

100 

100 

50 


The interpretation of these results is discussed for the case 
of a hypothetical sample of 100 random landmark pulses in the ,55ii 
channel data. In this case one would expect to find that 79 of the 
pulses suffered no degradation due to detector sharing (P(FD=1.0) = 

.79), and that the remaining 21 pulses suffered a degratation by a 
factor greater than 0 (1. -RPH (FD=1) and less or equal to 50 percent. 
Eighty three of the pulses, on the other hand, would be expected to 
have pulse amplitudes that were equal to or greater than 95 percent of 
the maximum. Consequently, one would expect to find that the number 
of pulses which suffered an amplitude degradation of greater than 5 
percent [100 (1-RPH (FD=. 9) ] was 17 or 100 (1-P (FD=. 9) ) . In a like 
manner one finds that the probability of a degradation or decrease 
greater than 25 percent is 8 percent (or 8 .pulses out of 100) , the 
probability of a decrease greater than 37 percent decrease is 4 per- 
cent, and the probability of a decrease greater than SO percent is zero. 
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A-2 SAMPLING DEGRADATION STATISTICS 


The discrete sampling of MSS imagery signals (as described in 
Section 2. 1.3. 3) can under certain conditions result in a decrease in 
the recorded pulse amplitude of a beacon landmark signature, correspond- 
ing to a fractional relative pulse amplitude (RPH) of less than unity. 
These conditions, and the statistical degradation that results when they 
exist are described as follows. 


With reference to the definitions and notations defined in the 
preceding discussion, we recall that during the time interval between 
successive samples of a given detector, the image spot translates through 
a distance AX=gL in the image plane, where g = T /DT is the ratio of , 
the sample repetition period (T^^) to the detector dwell time{DT). The 
data is said to be oversampled if T^^ <.DT, making g=AX/L<1.0. For 
the Landsat MSS, the value of g was shown (in Section 2, 1.3. 3) to be 
approximately equal to 0.74, The digital representation of an analog 
pulse signature will therefore generally consist of several non zero 
samples or pixels representing a particular distribution of intensities 
which depends on the exact but random value of the first non zero sam- 
ple and the relative values of AX, the detector width L, and the image 

spot diameter D , or, g e ax/L and a= D/L. 

s s 


■ The various possible pulse shapes associated with a given value 

of g and a are illustrated in Figure A-2 for g = 0.6 and for three 

values of a (1, 1/3, and 1/5) . The various pulse shapes for a given 

g and a were created by siit^sly varying the initial amount of the spot 

overlap (FDE X/D ) associated with the first non zero sample. This 
s 

variation of the initial condition was accomplished in discrete in- 
crements corresponding to = kL/lO (for k = 1, 2, 3. . .), and the' 
resulting (kth) pulse- shape was then obtained from the sequence of the 
(jth)” sampling positions defined by. 


Xkj " (0.6L)j (for j = 0,1,2,. . .), 


(A-5) 


which were used to determine the corresponding overlap 
and the corresponding RPHj, j ^^^kj ^ values from Equation 


fractions FD 


kj 


A-3. 
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Figure A-2 Efficiency and Shapes of Beacon Landmark Pulse Signatures 

for a Uniform Distribution of Equally Probable Initial Conditions 
Shown for 3 Fractional Spot Diameters (a = D /L) and a Sampling 
Step Size AX = 0.60L ® 






From Figure A-2 it is clear that the possibility of pulse ampli- 
tude degradation due to sampling is enhanced by large values of a and/ 
or g. The condition which guarantees an undegraded signal amplitude 
for at least one of the samples of a pulse signature is expressed by 
the relationship 

(g + a) <1.0. (A- 6) 

The validity of this condition is supported by the constructed pulse 
shapes shown in Figure A-2, for the two smaller values of a, which 
correspond to values of 4/5 and 28/30 for (g+a). Failure to satisfy 
this inequality condition does not mean that one of the RPH values can- 
not be unity, but only that a unity value is not guaranteed. This is 
illustrated for the third case, where a = 1.0, and it is seen that an 
undegraded pulse amplitude results for one point in the initial sample 
spectrum even though g+a is as large as 1.6. In this case the maximum 
pulse amplitude degradation, corresponding to another point in the 
initial sample spectrum, can be as large as 25 percent. Under these 
conditions one would expect the degradation of observed pulses to be 
randomly distributed over a range of from 0 to 25 percent of the maxi- 
mum possible pulse amplitude. 

From Figure A-2 it is clear that the maximum degradation (for the 
case when g+a <1.0) occurs for the initial sample condition that results 
in a symmetric pulse shape characterized by two equal relative-maximum 
sample signals, while the undegraded pulse shape (for a £ 1.0) is 

associated with a symmetric pulse shape having only one maximum sample 
signal. Thus, from the symmetry involved, one finds that the lowest 
relative pulse amplitude RPH* (FD*) for the maximally degraded pulse 
shape corresponds to a' fractional overlap of magnitude 

FD* = .5 [1 + (1-g) ] . (A-7) 

a 

For the general case then, in which g + a>1.0, one can use Equa- 
tions A-7 and A-3 to determine RPH* and the fractional degradation 
range DR* defined by DR* = (100-RPH*) percent of the maximum (under- 
graded) pulse amplitude. Thus, taking g = 0.74 and a^ = .21, .25, .29, 
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and .37 for the Landsat MSS, we conclude from the above agruments that 
the corresponding theoretical degradation ranges are 0, 0, 2, 

and 9 percent of the maximum possible pulse amplitudes in the corres- 
ponding wavelength channels. From this we conclude that the statisti- 
cal effects of pulse amplitude degradation due to discrete sampling 
(in the Landsat MSS) would not be appreciable compared to the nominal 
background noise associated with the ordinary scan signal. Thus in 
relation to the analysis of the results of the Landsat/Searchlight 
landmark experiments (discussed in Section 2.5), it appears that the 
, only meaningful conclusion which can be drawn from the above findings is 
that they serve only to magnify the conservative nature of the theor- 
etical pulse amplitude model which underestimated the observed near IR 
pulse amplitudes, on the average, by 10 percent. 
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APPENDIX B 


SOLAR ZENITH ANGLE VARIATIONS 
AT SUBSATELLITE POINTS AS A 
FUNCTION OF TIME AND LATITUDE 

Because EOS orbits are designed to be sun synchronous, the satellite will 
pass over the same point on the earth at essentially the same local time every 18 
days. However, even though the local time remains essentially the same, changes 
in the north/ south position of the sun with time of year (or season) will alter the 
solar zenith angle (0^) at the subsatellite point (as defined in Figure B-1), thus 
producing variations in the lighting conditions under which EOS imagery is obtained. 
The value of 0^ at any given subsatellite point can therefore be represented as a 
function of three parameters; the latitude of the subsatellite point, the local time at 
the descending node. (north to south equatorial crossing), and the time of year. For 
any given EOS orbit, however, the local time at the descending node remains 
essentially constant throughout the entire mission. 

A general expression for the variation of 0 with time of year and latitude 
is derived with the aid of Figures B-2 and B-3 as follows. Figure B-2 shows the 
sun and satellite position vectors S and R in an earth centered descending node 
reference frame (XYZ) which is regarded as fixed in inertial space. The x-axis 
defines the fixed direction from the center of the earth to the descending node, and 
the z-axis defines the fixed direction of the earth's spin axis. The satellite position 
vector (R) is referenced to the XYZ coordinate frame through the spherical polar 
coordinate angles 0 and (p, while the sun position vector (S) is referenced to the 
XYZ frame through the solar elevation angle (measured from the equatorial plane) 
and the descending node - sun angle y measured in the equatorial plane between the 
descending node and the projection (| ) of S onto the equatorial plane. The solar 
elevation angle (fj)) varies bi-yearly between plus and minus 23.5 degrees (in a 
manner to be described) and the descending node-sun angle (y) for an A.M. crossing 
time Tj^^ (measured in hours) is given by 

y = 15 (12-Tj^j^) degrees, (B-1) 

From Figure B-2, the components of the unit satellite and sun vectors 
(R and S) along the directions x, y, and z are clearly given by 

R = Sin (0) Cos (^) x + Sin (0) Sin (to) y + Cos (0) z, ■ (B-2) 

S = Cos (0) Cos (y) X 4 - Cos (01 Sin (y) y 4 - Sin (0) z. (B-3) 


B-1 



SATELLITE 




Figure B*1 . Solar Zenitii Angle Variations at Subsatellite Point 

as a Function of Season and Latitude for a 9:30 am 
(Descending Node) Sun Synchronous Orbit- 
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Figure B-2. Sun and Satellite Position Vectors in 

the Earth's Descending Node Reference Frame. 
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Figure B-3. Orientation of the Earth's Reference Frame (X,Y,Z), 
the Noon Equatorial Local Vertical ^ ,.and thp Solar 
Position Vector S in the Ecliptic Reference Frame 
(Xs-Ys,Zs). 
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Using the'vector relation'R •• S' =' cos we obtain a general expression for '0^ in 
the form: , . . ■ • 

0^ = arc cos [ cos (0) cos (y) cos (<p^ sin (0) 

+ cos ( 0')‘ sin'Cy)- sin {(p) sin (0) + sin (0) cos-('0)] (B-4) 


This equation as it stands is unsuitable for determining 0^ however because of the 
fact that the azimuth angle 0 is not an independent parameter. The angle (pcan be 
expressed, however^ in terms of the independent parameters 0 and the tilt angle . 
(6^y) of the orbital plane about the negative descending node axis (-X), measured 
relative to the Z-axis. This is accomplished by temporarily introducing another 
angle parameter (6 ) which defines the angle between the :satellite position vector R 
Sind the X-axis. Then from Figure B-2 it is clear that the components of R along Y 
and Z can also be expressed in terms of 0 • and 6, giving: 

zy 

' R = R ■ Y = Sin (0) Sin (0) = Sin'(6 ) Sin (6^ ) (B-5) 

y “ ^y 

and , . , 

= R ■ z = Cos (0> = Sin (6) Cos J (B-6) 

z zy 

Elimination of Sin (6 ) from the equations then gives the desired relation 
Sin (0) = tan'(0 ) / tan (9) (B-7) 

?y 


Defining j9 = 9 and Oi = 

zy 

for 0 in the form 
o 


ff /2-0 (for convenience) we obtain a general expression 


0 

o 


arc cos J cos (ij>) cos (y) cos (oi ) 


/l - tan^ (jg) \ 

V tan2 (9) / 


1/2 


tan-^ (9) 


+ Sin(y) Sin (o) tan (/5) 


+ Sin (0) Sin (a) j (B-8) 

which gives ® ^ in terms of the latitude (a), the solar elevation (0), and the fixed 
orbital parameters (yand j9). 


The explicit variation of 0with time is derived with the aid of Figure B-3 
which depicts a non-spinning earth coordinate frame (xyz) with the earth's equator- 
ial (xy) plane tilted through an angle 0 ^ = 23, 5 degrees in an ecliptic reference 

frame (X , Y , Z ) which is similarly regarded as being fixed or at rest in space, 
s s s 
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The angular position of the sun vector (S) relative to in the ecliptic plane 
at time (t) is then denoted by (t), with the vernal equinox being' defined for 
0 (t) = ff. The solar elevation angle i|i shown in Figure B-2 is here depicted as 

ij) (t) such that 0 (t) = -0, with C still denoting the perpendicular projection of S onto 
the equatorial plane (xy). The vector representing the noon equatorial local 
vertical lies in the plane-determined by Z and S, and the vectors {Xg>X,|^, and Y) 
all lie in the equatorial plane with the vector Y being in the Y^ plane and in- 
clined at an angle 0 ^ relative to Y^. 

The variation of 0 (t) with respect to the nearly uniform increase of 0 ^ (t> 
is derived from an examination of the direction cosines or components of a unit sun 
vector S in the XYZ frame. 

S = Cos (0 ) X + Sin (0_) Cos (\j) ) Y + Sin (0) Z (B-9) 

S s o 

From the normality condition ^ * S = 1 it follows that 

Cos^ (0 J + Sin^ (0 J Cos^ (0J. + Sin^ (0) = 1 . (B-101 

s so 

and Sin (0) = Sin (0.) Sin {0 K (B-H) 

Thus, 0 = - 0 (t) = - arc sin [ Sin (0 ) Sin (0 (t) ) ] defines the values of 0 

o s 

used in Equation B-8 to compute the yearly variation of 6 ^ shown in Figure B-1, 
for latitudes of 30 and 50 degrees north, a descending node crossing time 
.of 9:30 A.M. (y= 37. 5 degrees), and an orbital inclination of 99 degrees (/3 = 9 de- 
grees). 
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APPENDIX C 


DESIGN CONSIDERATIONS FOR 
SOLAR BEACON LANDMARKS 


In this appendix we discuss some of the theoretical and practical 
considerations relating to the use of curved mirrors as the effective 
radiation element in an artificial landmark station. From a theoreti- 
cal viewpoint we show that curved mirrors could be used to do the job, 
but that the use of such would require the costly addition of an auto- 
mated, two- axis, sun-tracking and mirror-pointing subsystem or, an un- 
commonly large curved mirror, in order to compensate for the large sea- 
sonal variations in the solar zenith angle (6^^) shown in Figure C-1. 

We note parenthetically, that neither of these alternatives appears 
to be as attractive, from the viewpoint of’ cost and reliability, as the 
alternative of using the compact and commercially available xenon search- 
light described in Section 2.6.3 of this report. In addition to provid- 
ing a theoretically adequate fixed beam radiation source for a nominal 
cost, a Searchlight Landmark Station (SLS) offers additional reliability 
and versatility through its capability to function at night and in every 
case where the overhead visibility is relatively good,' even though the 
slant range solar radiation may be highly attenuated due to haze and/or 
distant clouds. Also, because of the extremely small duty cycle required, 
the cost of operating a 4 kw searchlight for a period of just six min- 
utes every eighteen days can be considered insignificant in comparison 
to the overall operation costs for a national system of beacon land- 
mark stations. This is especially clear when one realizes that a real- 
istic design for the radiation element, and the landmark station as a 
whole, must include protection against damage that would otherwise re- 
sult from harsh environmental conditions. 

The large yearly variations in the sun's zenith and azimuth posi- 
tions are therefore shown 'to introduce technical complications which 
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Figure C-1 Solar Zenith Angle Variations at Subsatellite Point 

as a Function of Season and Latitude for a 9:30 am 
(Descending Node) Sun Synchronous Orbit* 
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undermine the seemingly ideal virtures of using mirrors for the effect- 
ive radiation element in a national system of beacon landmarks. These 
complications are not such as to make a Mirror Landmark Station (MLS) 
impractical but only less desirable (theoretically) than an SLS from 
the viewpoint of reliability, cost, and versatility. Although the ex- 
periments described in Sections 2.4 and 2.5 have confirmed the anti- 
cipated theoretical attributes of an SLS , it is nevertheless of inter- 
est to examine the theoretical performance and design considerations 
relating to the reflector component of an MLS. 

In the following analyses we assume that the satellite passes 
directly over the landmark and we estimate the performance and design 
characteristics of both relatively small curved mirrors , which would 
require a sun-tracking and mirror pointing subsystem, and relatively 
large curved mirrors which could be maintained in a fixed orientation 
relative to the earth. Estimates of the actual mirror sizes in each 
case are determined on the basis of the intensity criteria described 
next., 


C.l , SOLAR BEACON INTENSITY CRITERIA 
* 

The -basic intensity criterion imposed on any landmark beacon is 
that it be capable of generating a unique and readily detectable pulse 
signature in the EOS MSS imagery data, when the interception is within 
a 1.5 by 6.0 degree cone. Automated recognition of a solar beacon land- 
mark, as has been demonstrated with xenon searchlight landmark data, 
can be achieved by using a combination of \miqueness criteria based on 
the known spectral shape of the beacon landmark radiation and a set of 
spectral pulse amplitude thresholds for the MSS channels. 

In the case of .a xenon SLS for example, the xenon flair (shown 
in Figure C-2) , which is 'spectrally confined within the near IR chan- 
nel of the Landsat and EOS MSSs, provided a natural spectral uniqueness 
criterion which allowed low amplitude xenon searchlight pulses to be 
readily detected within an otherwise noisey background of the numerous 
and relatively large solar beacon pulses that exist in the MSS imagery 
data, due to the variety of natural and man made specular reflectors 
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Figure c-2 Normalized Spectral Distributions for Solar 
Radiation {Upper Atmosphere) and Xenon 
Searchlight (in the kilowatt range). * 



that exist on the earth's surface. Under these conditions it was possi- 
ble to construct an automated SLS detection scheme with a high degree 
of selectivity (or noise rejection) for a pulse-amplitude threshold 
sensitivity of only 7. percent (of full scale) ’in the near IR channel 
of the Lands at MSS. 

A practical searchlight or solar beacon landmark should 'be capable 
of generating a pulse in the most sensitive or dominant MSS channel with 
an amplitude that would, under normal conditions, be much larger than 
the threshold sensitivity of the detection scheme. The safety factor 
must be large enough to allow for a reduction in the recorded pulse 
amplitudes by a factor of two, due to the phenomenon of detector shar- 
ing, as well as likely reductions in the beam intensity due to ageing 
and/or contamination of the optics, and also for possible increases in 
the optical thickness .of the atmosphere relative to the clear standard 
atmospheric optical thicknesses assumed in the theoretical model. In 
addition, an allowance has to be made for the presently unknown pulse- 
an^litude threshold sensitivity of a near real time, on board, analog 
beacon landmark detection filter. This analog threshold sensitivity 
could conceivably be two to three times larger than the pulse amplitude 
threshold sensitivity for a digital filter due to the fact that the 
uniqueness test has to be performed in near real time. 

Based on, these considerations , a prototype xenon searchlight for 
an EOS SLS was selected with the capacity to generate a nominal (un- 
degraded) pulse amplitude of 46 percent of full scale, which is over 
six times larger than the demonstrated pulse amplitude sensitivity of 
the automated digital pulse amplitude detection scheme discussed in 
Section 1.6. 

The corresponding design goal for the nominal (undegraded) ampli- 
tude of a solar beacon landmark pulse cannot be accurately estimated 
without first having some idea of the corresponding minimum threshold 
sensitivity that could be achieved in a highly selective automated 
detection scheme. A high degree of discrimination against other 
naturally occuring solar beacon signatures could conceivably be achiev- 
ed by employing a transmission filter to limit the spectrum of the 
reflected light from an MLS to the range associated with just one 
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specific MSS channel. It will here be assumed that the use of such a 
filter would provide an adequate MLS detection selectivity for a pulse 
amplitude threshold of 7 percent of full scale in the .75 micron chan- 
nel of the EOS MSS. Thus with due allowance for the attenuation of 
the sunlight due to the filter, the design goal for solar beacon 
landmark signature will be assumed to be the same as that for a xenon 
searchlight landmark signature. Namely that the nominal (undegraded) 
pulse amplitude be 46 percent of the full scale range for the .75 micron 
channel of the EOS MSS, over an interception cone angle of 1.5 by 6.0 
degrees . 

C.2 THEORETICAL SIGNATURE CHARACTERISTICS OF SOLAR BEACON LANDMARKS 

In this section we develop a theoretical model for predicting the 
amplitude of a solar beacon landmark pulse in terms of the size of a 
curved mirror and the angular position of the sun relative to the local 
vertical. The model, subject to the aforementioned intensity criteria, 
is then used to determine the approximate size of a standardized mirror 
for an EOS MLS. 

C.2.1 General Geometric Considerations 


The function of an MLS mirror, as illustrated in Figure C-3 , is 
to reflect incident solar radiation into a small cone which can be 
maintained in a fixed orientation relative to the earths surface, 
irrespective of the seasonal changes in the sun’s position. The orient- 
ation being that which causes the center of the cone to intercept the 
nominal (time-average) orbital path of the satellite at right angles. 

In addition, the cone size or beam spread must be sufficient to allow 
for the maximum anticipated deviations in the sighting angle due to 
deviations of the spacecraft attitude and orbital position from nomi- 
nal. In the present case it has been assumed that this requirement 
can be met by a beam spread of 1.5 degrees in the intrack direction 
and G.O degrees in the crosstrack direction. 
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These conditions could be satisfied in principle by employing a 
fixed curved mirror, as illustrated in Figure C-4 which is large 
enough to accommodate the seasonal changes in the sun's position for 
all latitudes in the jj.S. and for the time of day determined by the 
particular sun syndronous orbit employed. The seasonal variations of 
the sun's zenith angle throughout the U-.S., as shown in Figure C-1 
(for the case of a 9:30 A.M. decending node), suggest that the angular 
dimensions of a standardized stationary mirror would have to be large 
enough to accommodate a 50 degree range of solar zenith angles, with 
a corresponding azimuth angle range of about 54 degrees. The physical 
size of such a mirror would then be determined by the magnitude of the 
radius of curvature that is required to satisfy the intensity criteria 
for an MLS. 

These considerations illuminate an important physical character- 
istic of a spherical convex mirror. Its three-dimensional size or 
volume is directly proportional to its surface area ) which is 

fixed or determined by the product of two functionally independent 
design parameters: the solid angle subtended which is determined 

entirely by the angular zenith and azimuth ranges of the incident 
radiation that must be accommodated, and the radius of curvature (r) 
which is determined entirely by the reflected intensity requirements 
for a given illumination intensity. This characteristic of a curved 
mirror tends to conflict with the intuitive notion that the -reflected 
intensity should depend explicitly on area of the mirror surface, which 
is only true for the special case of a flat mirror. 

Thus, if the size of a standardized stationary MLS mirror would 
be impractically large, as the present analysis will suggest, then the 
only ■ alternative left is to minimize the ranges of the solar zenith and 
azimuth angles that have to be accommodated through the use of an 
automated two-axis sun tracking and mirror pointing control system. 

Thus minimizing the values of the mirror's polar and azimuth angle in- 
crements A0 and A<j),and hence its overall size for a given radius of 
curvature. The determination of these minimum polar and azimuth angles 
is discussed next. 
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Figure C-4 Geometry and Size Parameters for a 
Standardized Stationary Curved Mirror 
Landmark Element. 
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C.2.2 Minimum Polar and Azimuth Angle Ranges for a Curved Solar 
Beacon Mirror 


The minimum values of the curved mirror’s solar and azimuth angle 
increments, A0 and A(j>, required to provide a reflected beam spread of 
1.5. by 6.-0 degrees along the iritrack and cross track directions, res- 
pectively, were determined with the aide of two CSDL computer programs. 
One program (REV.ALIS) was used to compute the seasonal changes in the 
zenith and azimuth angles of the sun, relative to both a longitude- 
latitude surface coordinate frame and an intrack-crosstrack surface 
coordinate frame for latitude points between 30 and 50 degrees north. 
Then 'using the crosstrack direction as the azimuth reference, it was 
found that the solar zenith and azimuth angles varied over the ranges 
of 30 to 80 degrees and +2 to -50 degrees, respectively. The second 
computer program (REV. SPHERE) was then used to compute the minimum 
values of A0 and Ac}) required to assure a minimum intrack-crosstrack 
beam spread of 1.5 by 6.0 degrees over these solar zenith and azimuth 
ranges for the case where a gimbaled two-axis, sun-tracking and mirror- 
pointing system is employed to track the sun’s position and to keep 
the reflected beam oriented in the vertical direction. The results 
showed that the angles A0 and A4> subtended by the mirror would have to 
be on the order of 4.6 and 2.0 degrees respectively, for a nearly rect- 
angular mirror segment, corresponding to a mean polar angle of 6^(MAX)/2 
= 40 degrees. 

Thus, for a given radius of curvature, these curved mirror design 
parameters serve to determine the overall size of the curved mirror 
element for a standardized MLS that could be located at any latitude 
within the continental U.S., and at any point in a crosstrack strip 
which falls within about +10 nautical miles of the nominal subsatellite 
point ground track. These figures are based on an assumed sun tracking 
and mirror pointing accuracy of approximately +0.1 degrees in the res- 
pective zenith and azimuth angles. The minimum radius of curvature 
required to assure a nominal (undegraded) MLS pulse amplitude of 46 
percent in an EOS MSS is discussed next. 
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C.2.3 The Minimxm Radius. of Curvature for a Curved Solar Beacon Mirror 


The nominal (undegraded) beacon pulse amplitude (BPA^) in the 
ith channel of an MSS corresponding to a solar beacon can be. .computed 
as a fractional percent of the full scale response for the ith channel 
by means of the equation 


BPA^ = 100 PB^/PEM^ %, 


(C-1) 


where PB^^ is the solar beacon power contained in the ith channel image 
spot (or blurr circle) in the MSS image plane (for further details on 
this see discussion of Beacon Landmark Imaging Characteristics in 
Section 2.1.3), and PEM^ is the maximum recordable extended source 
power impinging on a given detector in the ith channel. PEM^ is there- 
fore the detector power associated with a maximum (100 of percent full 
scale) scan signal in the ith channel of the MSS. In terms of the MSS 
parameters 


PEDf 


^i 


gmi D A 
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where k^ is the MSS optical transmission factor for the ith channel; 

BR . is the calibrated (ith channel) Bolometric Intensity Radiance 
gmi 

[watts’/M— str] at the MSS aperture associated with a maximum or full 
scale scan signal in the ith channel; ^2^=9^ is the solid IFOV of the 
MSS, given ‘by the square of the linear IFOV (0j^) ; and A^ is the 
aperture area. (For further details on this see discussion of Extended 
Source Imaging Characteristics in Section 2.1.2). 


The solar beacon image power on the other hand can be expressed 


as 


PB. = tjPE^ 3W. a^, 


(C-3) 


where PR^ is the Power Radiance [watts/v^^-str] at the aperture associa- 
ted with the solar beacon (averaged over the ith bandwidth); BW. [y] is 

2 . 

the bandwidth of the ith channel; and Q^=A^/R is the solid angle sub- 
tended by the MSS aperture at an altitude R. The solar beacon power 
radiance is further defined by the expression 

-OT 



where PR? is defined as the unattenuated power radiance of the. solar 
beacon; OT^ is the vertical optical thickness of the atmosphere (aver- 
aged over the ith bandwidth) ; 0^ is the solar zenith angle as defined 
in Figure C-1, and is the MLS fil-ter transmission factor. The 
terra exp [-OT^ (1+sec (0^) ) ] represents the combined vertical and slant 
range atmospheric transmission factor based on a flat atmosphere approx- 
imation. 

The unattenuated power radiance of solar radiation reflected from 
a curved mirror is derived with reference to the ray geometry shown in 
Figure C-5, Radiation propagating in the -Z direction, and assvraied 
initially to be plane parallel, is intercepted by a small spherical 
reflecting surface of area whith subtends a small solid angle 

about a polar angle 6 = If Id) is the irradiance or 

spectral intensity [watts/m~p] of the incident radiation, then the 
spectral power density [watts/jj] of the vertically reflected beam is 
given by 


PD(X) = Id) Ajjjcos(0), (C-5) 

= r^ I(X)f2jjj cos(0). 

A ray which impinges on the mirror surface at a point P(r,0,({i) will 
have an angle of incidence equal to 0 and will be reflected through an 
angle of 20 in the azimuth plane (j) = constant. The polar angle of a 
reflected ray (Qq) is therefore related to the polar angle 0 of its 
point of origin on the mirror surface by the transformation 0j^=20. 

The solid angle of reflected beam 

= Sin(0j^)A0^ Aij) , (C-6) 

is therefore related to the solid angle of the mirror 

= Sin (0)ABA4) , (C-7) 

by = Sin (20) (2A0)A0 (C-8) 

= 4 cos (0) Sin (0 ) A6At() 

= 4 cos(0)Q • (for << 1.0) 

M M 
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If the incident radiation is not exactly parallel however, but 

is derived from an extended source like the sun, then the solid angle 

of the reflected radiation will be increased by the solid angle 

=ire^/4 defining the solid divergence of a uniformly bright disc 
s s 

which subtends an angle 6^ = .53 degrees relative to a point on the 
earth's surface. Thus 


£2 = 4 cos ( 0 ) 


(C-9) 


and the unattenuated power radiance (PR?) of solar radiation reflected 
from a curved mirror can then be epxressed in the form 


PR? (9) = ^^i 


r'^Ii £2 j4 Cos (0) 
4cos (0 ) 


(C-10) 


which can be simplified to 


PR? = PR?(0) 




4 (1+S) 


, watts ^ 

\ 1 r 

y-str 
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where 


£2 


4 cos { 0 ) £2 


m 
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PR? is therefore expressed rather simply in terms of the radius of 
curvature (r) of the mirror; the solar irradiance at the top of the 
atmosphere (averaged over the bandwidth of the ith MSS channel) ; and 
the parameter 5 . Since 9 is restricted to the value 0^/2 in order to 
provide a vertically oriented beam, 6 can be expressed in the equi- 
valent form 


6 


" ^s 

8 Sin (6 q) A9A(j) 


(C-13) 


Thus for the conditions stated in Section C.2.2, (0^=0. 53 degrees, 

9^{MAX)=80 degrees, A0=4.6 degrees, and A<1)=2,0 degrees) we have 
6(MIN)=.012 and 6 (MAX) =6 (9^=30 degrees)=.024. This shows that the 
changes in the angular orientation (6) of the mirror, as it tracks the 
sun’s position to maintain 9=0^/2, will have very little affect on the 
intensity of the reflected beam. 
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Finally, using Equations C-1 to C-4 and C-11, we obtain the following 
simplified' expression for the theoretical pulse amplitude of solar 
beacon landmark 


where 


with 


BPAi 



„ _ 25 I. BW- FT. AT. 

— X X J- X 

®^GMi 

AT. = iOTi[l+Sec(0^)] 


% Full Scale (C-14J 

(C-15) 


defines the MSS channel sensitivity constants. The values of shown 
in Table C-1 were computed for a satellite altitude (R) of 700 kilor 
meters, an MSS IPOV (Ojj) of 4 3 micro radians, a filter transmission , 
factor (FTj^) of 0.8 in all channels, and a solar zenith angle of 80 
degrees (5=.012), with the values of taken from Figure C-2. 


TABLE C-1 

EOS MSS SENSITIVITY CONSTANTS 
FOR A CURVED . MIRROR SOLAR BEACON 


EOS MSS CHANNEL 

1 

2 

3 

4 

Center Wavelength [p] ' 

.55 

.65 

.75 

.95 


.1 

.1 

.1 

. 3 

OT^E-] 

.379 

.295 

.242 

.195 

ATfE-] 

.077 

.136 

.195 

.268 

[watts/M^-y] 

1759 ■ 

1583 

1251 

812 

BRcjji [watts/M^-str] 

24.8 

20.0 

17.6 

46.0 

G^[%FS/M^] 

.01191 

.02348 

.03024 

.03095 

r^[M] 

62 

44 

39 

39 


The radius of curvature required to produce a nominal (undegraded) 
pulse amplitude of 46 percent of full scale was then obtained from the 
equation r^^ = (46/Gj^) ^/2 = 62, 44, 39, and 39 meters respectively. The 
corresponding minimum size of a standardized sun tracking MLS mirror 
would therefore be characterized by rectangular curvalinear dimensions 
on the order of 
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= r Sin(40)A4> = .875 M or 2.9 feet, . (C-16) 

and £q = rA6 = 3.13Morl0.3 feet. 

These mirror dimensions are larger than would be required to satisfy 
the design goal minimum pulse amplitude criterion (of 46 percent of 
full scale) at more elevated sun positions and/or at latitudes below 
50 degrees north. Consequently we emphasize that the above curved 
mirror dimensions relate specifically to a standardized EOS MIS mirror 
element design which theoretically satisfies the following design re- 
quirements : It permits MLSs to be located at any latitude within the 
continental U.S.A.,and with a sun-tracking and mirror-pointing sub- 
system, having an absolute accuracy of about +0.1 degrees in azimuth 
and elevation, the mirror element can maintain an intrack-cross track 
beam spread of at least 1.5 by 6.0 degrees {relative to the nominal 
ground' track) throughout the entire year. In addition, it provides a 
nominal (undegraded) pulse amplitude (for an EOS MSS IFOV of 43 micro- 
radians or less) under the most unfavorable conditions of sun position 
and/or latitude. 

The necessity of employing a sun-tracking and mirror-pointing 
device could of course be circumvented by utilizing fixed curved mirrors 
that are large enough to accommodate the seasonal changes in the sun's 
position at the time of flyby. With reference to Figure C-4 again, the 
size and approximate shape of such mirrors, consistent with a standard- 
ized design, would then be governed by the following values for the 
relevant angular parameters ; 


r 

= 128 feet 

0 (center) 

= 27.5 degrees 

A0 

— 28.0 degrees 

A'lp 

= 54.0 degrees 

h 

= 13.5 degrees 

02 

= 41; 5 degrees 
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The corresponding physical dimensions would then have the values : 

c = 63 feet 

W2 = 80 feet 

W1 = 28 feet 

a =55 feet 
h =29 feet 

Fixed mirrors, it seems, would not be desirable for an MLS design, even- 
if it could be shown that the overall cost per station would be less 
than that for a station design incorporating an automated sun-tracking 
and mirror-pointing device. In the future it may be possible to con-r 
sider smaller and more practical mirror systems based on a different 
set of design requirements for such things as the minimum acceptable 
pulse amplitude, and/or the uncertainties in the satellite's attitude 
and orbital position, and/or the IFOV of future EOSs, Any improve- 
ments in this regard would of course serve to reduce the size and cost 
of the sun-tracking and mirror-pointing subsystem required for an MLS 
and, also, the size and cost of the searchlight component required for 
an SLS. 
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